Packet 7 (KW first pass)

Last updated: 2019-04-15

This is Kara’s first pass at looking at the Packet 7 data (with data from 2019-03-18). (Note: PDF printed portrait, 80% zoom, minimum margins.)

Confirmatory hypotheses (as specified in preregistration)

NOTE: I included all participants, even if they failed attention checks, since we did not specify that we would exclude participants in the preregistration. At some point we should try running everything on the subset of participants who passed all attention checks just so that we know, but this decreases our sample sizes pretty dramatically.

ALSO NOTE: I included only the second sample of Ghanaian participants, who were recruited according to our standards and as specified in the preregistration. For the most part, all results are the same when including the full, “double” sample of Ghanaian participants.

Hypothesis #1

Hypothesis #1: People who score higher on the Tellegen Absorption Scale (Tellegen & Atkinson, 1974) also report more spiritual experiences.

Preregistered analysis: “We will standardize scores on each measure, collapsing across field sites, and then conduct a mixed-effects linear regression, using absorption scores to predict reports of spiritual experiences, allowing for differences across our five field sites and across the two scales we used to measure spiritual experience. Using the “lme4” package for R (Bates, 2018), we will specify this model as follows: experience ~ absorption * site * scale + (1 | subject), where “scale” is a categorical variable comparing the two spiritual experience scales (effect-coded), and “site” is a categorical variable comparing the five field sites (effect-coded). Our prediction is that the main effect of absorption on experience will be significantly positive (b > 0, t > 2, p < 0.05 using Satterthwaite’s degrees of freedom method as implemented in the “lmerTest” package for R, Kuznetsova, Brockhoff, & Christensen, 2017). (The rest of the parameters included in the model will be of great interest to us, but are exploratory.)

If we have trouble fitting these models or there are other signs of model fit problems, we will consider the following alternative model: experience ~ absorption * scale + (1 | site/subject)."

Here are our results:

parameter	b	standard error	df	t	p	significant
(Intercept)	0.03	0.03	495	1.12	0.263
abs	0.21	0.03	495	6.55	<0.001	*
p7_ctry1	-0.32	0.06	495	-5.44	<0.001	*
p7_ctry2	0.79	0.06	495	13.58	<0.001	*
p7_ctry3	-0.11	0.06	495	-1.78	0.075
p7_ctry4	-0.95	0.06	495	-15.26	<0.001	*
scale1	0.00	0.02	495	-0.17	0.863
abs:p7_ctry1	0.04	0.05	495	0.72	0.474
abs:p7_ctry2	-0.01	0.07	495	-0.18	0.86
abs:p7_ctry3	0.12	0.07	495	1.79	0.073
abs:p7_ctry4	-0.04	0.07	495	-0.60	0.551
abs:scale1	-0.03	0.02	495	-1.52	0.129
p7_ctry1:scale1	-0.06	0.03	495	-1.86	0.064
p7_ctry2:scale1	0.04	0.03	495	1.14	0.255
p7_ctry3:scale1	0.07	0.04	495	1.99	0.047	*
p7_ctry4:scale1	-0.18	0.04	495	-4.90	<0.001	*
abs:p7_ctry1:scale1	-0.03	0.03	495	-0.88	0.381
abs:p7_ctry2:scale1	-0.03	0.04	495	-0.68	0.499
abs:p7_ctry3:scale1	0.00	0.04	495	-0.01	0.992
abs:p7_ctry4:scale1	0.04	0.04	495	0.87	0.387

We found evidence in favor of Hypothesis #1: As predicted, the main effect of absorption was significantly positive (see bolded row in table above).

For every increase of 1 standard deviation in absorption, participants (on average) increased by 0.21 standard deviations in their spiritual experiences. In concrete terms, this means that for every increase of 6.51 points on the Absorption scale (range: 0-34), participants (on average) increased by 3.46 points on the Daily Spiritual Experiences scale (range: 0-70) or 3.16 points on the Spiritual Events scale (range: 0-80).

Hypothesis #2

Hypothesis #2: People who score higher on our two porosity indices also report more spiritual experiences.

Preregistered analysis: “We will standardize scores on each measure, collapsing across field sites, and then conduct a mixed-effects linear regression, using porosity scores (from each of our two indices) to predict reports of spiritual experiences across the 2 spiritual experience scales included in our survey packet, allowing for differences across our five field sites and across the two scales we used to measure spiritual experience. Using the “lme4” package for R (Bates, 2018), we will specify this model as follows: experience ~ porosity * site * scale + (1 | subject), where “scale” is a categorical variable comparing the two spiritual experience scales (effect-coded), and “site” is a categorical variable comparing the five field sites (effect-coded). We will run two versions of this model, one using each of the two “porosity” measures described above. For each regression, our prediction is that the main effect of porosity on experience will be significantly positive (b > 0, t > 2, p < 0.05 using Satterthwaite’s degrees of freedom method as implemented in the “lmerTest” package for R, Kuznetsova, Brockhoff, & Christensen, 2017). (The rest of the parameters included in the model will be of great interest to us, but are exploratory.)

If we have trouble fitting these models or there are other signs of model fit problems, we will consider the following alternative model: experience ~ porosity * scale + (1 | site/subject)."

Here are our results, using the Porosity scale (Dulin):

parameter	b	standard error	df	t	p	significant
(Intercept)	0.07	0.04	495	1.80	0.073
por	0.38	0.04	495	10.20	<0.001	*
p7_ctry1	-0.19	0.06	495	-3.17	0.002	*
p7_ctry2	0.58	0.07	495	8.56	<0.001	*
p7_ctry3	-0.17	0.06	495	-2.74	0.006	*
p7_ctry4	-0.59	0.08	495	-7.34	<0.001	*
scale1	0.01	0.02	495	0.39	0.694
por:p7_ctry1	0.20	0.06	495	3.25	0.001	*
por:p7_ctry2	-0.08	0.07	495	-1.17	0.245
por:p7_ctry3	0.00	0.07	495	0.05	0.963
por:p7_ctry4	0.04	0.08	495	0.54	0.59
por:scale1	0.00	0.02	495	0.14	0.89
p7_ctry1:scale1	-0.08	0.04	495	-2.06	0.04	*
p7_ctry2:scale1	0.11	0.04	495	2.48	0.014	*
p7_ctry3:scale1	0.07	0.04	495	1.88	0.06
p7_ctry4:scale1	-0.12	0.05	495	-2.40	0.017	*
por:p7_ctry1:scale1	-0.03	0.04	495	-0.75	0.456
por:p7_ctry2:scale1	-0.14	0.04	495	-3.29	0.001	*
por:p7_ctry3:scale1	-0.01	0.05	495	-0.15	0.878
por:p7_ctry4:scale1	0.08	0.05	495	1.62	0.106

And here are our results, using the Martha Story scale:

parameter	b	standard error	df	t	p	significant
(Intercept)	0.06	0.03	495	1.94	0.053
mm	0.26	0.03	495	7.60	<0.001	*
p7_ctry1	-0.18	0.06	495	-2.95	0.003	*
p7_ctry2	0.66	0.06	495	10.31	<0.001	*
p7_ctry3	-0.26	0.06	495	-4.58	<0.001	*
p7_ctry4	-0.75	0.07	495	-10.67	<0.001	*
scale1	0.01	0.02	495	0.39	0.693
mm:p7_ctry1	0.20	0.06	495	3.57	<0.001	*
mm:p7_ctry2	-0.07	0.06	495	-1.20	0.23
mm:p7_ctry3	-0.01	0.06	495	-0.14	0.888
mm:p7_ctry4	0.05	0.08	495	0.64	0.521
mm:scale1	-0.04	0.02	495	-2.14	0.033	*
p7_ctry1:scale1	-0.11	0.04	495	-3.11	0.002	*
p7_ctry2:scale1	0.10	0.04	495	2.66	0.008	*
p7_ctry3:scale1	0.07	0.03	495	2.10	0.036	*
p7_ctry4:scale1	-0.17	0.04	495	-3.96	<0.001	*
mm:p7_ctry1:scale1	-0.06	0.03	495	-1.83	0.068
mm:p7_ctry2:scale1	-0.09	0.04	495	-2.48	0.014	*
mm:p7_ctry3:scale1	0.01	0.04	495	0.28	0.778
mm:p7_ctry4:scale1	0.08	0.05	495	1.60	0.111

In both cases, we found evidence in favor of Hypothesis #2: As predicted, the main effect of porosity was significantly positive in both regression analyses (see bolded rows in tables above).

For every increase of 1 standard deviation in porosity as indexed by the Porosity (Dulin) scale, participants (on average) increased by 0.38 standard deviations in their spiritual experiences. In concrete terms, this means that for every increase of 7.47 points on the Porosity scale (range: 0-34), participants (on average) increased by 6.33 points on the Daily Spiritual Experiences scale (range: 0-70) or 5.79 points on the Spiritual Events scale (range: 0-80).

For every increase of 1 standard deviation in porosity as indexed by the Martha Story scale, participants (on average) increased by 0.26 standard deviations in their spiritual experiences. In concrete terms, this means that for every increase of 10.55 points on the Martha Story scale (range: 0-54), participants (on average) increased by 4.36 points on the Daily Spiritual Experiences scale (range: 0-70) or 3.99 points on the Spiritual Events scale (range: 0-80).

Hypothesis #3

Hypothesis #3: The positive relationship between absorption and experience is statistically significant in each field site individually.

Preregistered analysis: “For each of our five field sites, we will standardize scores on each measure (considering only participants from that field site), and then conduct the same analyses specified above for Hypothesis #1 to assess the relationship[] between absorption and spiritual experience […] in each field site individually. For Hypothesis #3, our models will take the following form: experience ~ absorption * scale + (1 | subject). […] Again, our prediction is that the main effect of absorption on experience and of each porosity measure on experience will be significantly positive (b > 0, p < 0.05).”

Here are our results:

parameter	b	standard error	df	t	p	significant
US
(Intercept)	0.00	0.09	100	0.00	1
abs	0.30	0.09	100	3.40	<0.001	*
scale1	0.00	0.04	100	0.00	1
abs:scale1	-0.06	0.04	100	-1.67	0.097
Ghana
(Intercept)	0.00	0.08	98	0.00	1
abs	0.20	0.08	98	2.42	0.017	*
scale1	0.00	0.06	98	0.00	1
abs:scale1	0.00	0.06	98	-0.06	0.955
Thailand
(Intercept)	0.00	0.07	109	0.00	1
abs	0.40	0.08	109	5.31	<0.001	*
scale1	0.00	0.04	109	0.00	1
abs:scale1	-0.03	0.05	109	-0.61	0.546
China
(Intercept)	0.00	0.09	98	0.00	1
abs	0.26	0.09	98	3.00	0.003	*
scale1	0.00	0.04	98	0.00	1
abs:scale1	-0.02	0.04	98	-0.41	0.684
Vanuatu
(Intercept)	0.00	0.08	90	0.00	1
abs	0.15	0.08	90	1.74	0.085
scale1	0.00	0.06	90	0.00	1
abs:scale1	0.01	0.06	90	0.14	0.892

We found evidence in favor of Hypothesis #3 in all sites except Vanuatu: As predicted, the main effect of absorption was significantly positive in the US, Ghana, Thailand, and China—but this was not the case in Vanuatu (although some would consider the effect “marginally significant”; see bolded rows in table above).

(See other relevant plots under “Hypothesis #1,” above.)

Hypothesis #4

Hypothesis #4: The positive relationship between absorption and experience is statistically significant in each field site individually.

Preregistered analysis: “For each of our five field sites, we will standardize scores on each measure (considering only participants from that field site), and then conduct the same analyses specified above for Hypothesis #2 to assess the relationship[…] between porosity and experience in each field site individually. […] For Hypothesis #4, our models will take the following form: experience ~ porosity * scale + (1 | subject). […] Again, our prediction is that the main effect of absorption on experience and of each porosity measure on experience will be significantly positive (b > 0, p < 0.05).”

Here are our results, using the Porosity scale (Dulin):

parameter	b	standard error	df	t	p	significant
US
(Intercept)	0.00	0.07	100	0.00	1
por	0.54	0.07	100	7.30	<0.001	*
scale1	0.00	0.04	100	0.00	1
por:scale1	-0.02	0.04	100	-0.49	0.628
Ghana
(Intercept)	0.00	0.08	98	0.00	1
por	0.29	0.08	98	3.64	<0.001	*
scale1	0.00	0.06	98	0.00	1
por:scale1	-0.06	0.06	98	-1.04	0.303
Thailand
(Intercept)	0.00	0.07	109	0.00	1
por	0.41	0.07	109	5.52	<0.001	*
scale1	0.00	0.05	109	0.00	1
por:scale1	0.00	0.05	109	0.10	0.917
China
(Intercept)	0.00	0.07	98	0.00	1
por	0.54	0.07	98	7.34	<0.001	*
scale1	0.00	0.04	98	0.00	1
por:scale1	0.05	0.04	98	1.14	0.257
Vanuatu
(Intercept)	0.00	0.08	90	0.00	1
por	0.19	0.08	90	2.25	0.027	*
scale1	0.00	0.06	90	0.00	1
por:scale1	0.11	0.06	90	1.83	0.07

Using the Porosity scale, We found evidence in favor of Hypothesis #4 in all sites: As predicted, the main effect of porosity was significantly positive in the US, Ghana, Thailand, China, and Vanuatu (see bolded rows in table above).

And here are our results, using the Martha Story scale:

parameter	b	standard error	df	t	p	significant
US
(Intercept)	0.00	0.08	100	0.00	1
mm	0.49	0.08	100	6.36	<0.001	*
scale1	0.00	0.04	100	0.00	1
mm:scale1	-0.11	0.04	100	-2.93	0.004	*
Ghana
(Intercept)	0.00	0.08	98	0.00	1
mm	0.20	0.08	98	2.42	0.017	*
scale1	0.00	0.06	98	0.00	1
mm:scale1	-0.10	0.06	98	-1.88	0.062
Thailand
(Intercept)	0.00	0.08	109	0.00	1
mm	0.33	0.08	109	4.17	<0.001	*
scale1	0.00	0.04	109	0.00	1
mm:scale1	-0.04	0.05	109	-0.83	0.411
China
(Intercept)	0.00	0.08	98	0.00	1
mm	0.38	0.08	98	4.68	<0.001	*
scale1	0.00	0.04	98	0.00	1
mm:scale1	0.00	0.04	98	0.04	0.968
Vanuatu
(Intercept)	0.00	0.09	90	0.00	1
mm	0.09	0.09	90	1.07	0.288
scale1	0.00	0.06	90	0.00	1
mm:scale1	0.03	0.06	90	0.48	0.632

Using the Martha Story scale, We found evidence in favor of Hypothesis #4 in all sites except Vanuatu: As predicted, the main effect of Martha Story was significantly positive in the US, Ghana, Thailand, and China—but this was not the case in Vanuatu (see bolded rows in table above).

(See other relevant plots under “Hypothesis #2,” above.)

Secondary exploratory analyses

“Secular” experiences

From preregistration: “First, in order to explore the question of whether affirmation of spiritual experience, porosity and absorption are simply due to the religiosity of participants, we will conduct additional analyses focusing on the two “secular” experience scales (the revised Launay Slade scale, Morrison, Wells and Nothard 2000]; and the Sheep-Goat [paranormal] scale, Thalborne & Delin, 1993). These analyses will be identical to the primary analyses outlined above, except that they will only include these two “secular experience” scales and not the two spiritual experience scales. We consider these to be exploratory analyses.”

Parallel to Hypothesis #1

parameter	b	standard error	df	t	p	significant
(Intercept)	0.01	0.03	495	0.34	0.731
abs	0.36	0.03	495	10.70	<0.001	*
p7_ctry1	0.14	0.06	495	2.26	0.025	*
p7_ctry2	0.24	0.06	495	3.74	<0.001	*
p7_ctry3	-0.05	0.07	495	-0.74	0.457
p7_ctry4	-0.23	0.07	495	-3.44	<0.001	*
scale1	0.00	0.03	495	0.06	0.949
abs:p7_ctry1	0.10	0.06	495	1.83	0.068
abs:p7_ctry2	-0.05	0.07	495	-0.75	0.451
abs:p7_ctry3	0.06	0.07	495	0.79	0.431
abs:p7_ctry4	0.04	0.07	495	0.50	0.617
abs:scale1	-0.01	0.03	495	-0.36	0.722
p7_ctry1:scale1	-0.25	0.05	495	-4.79	<0.001	*
p7_ctry2:scale1	-0.09	0.05	495	-1.62	0.107
p7_ctry3:scale1	0.19	0.06	495	3.37	<0.001	*
p7_ctry4:scale1	0.08	0.06	495	1.41	0.16
abs:p7_ctry1:scale1	0.00	0.05	495	0.00	0.999
abs:p7_ctry2:scale1	-0.05	0.06	495	-0.85	0.398
abs:p7_ctry3:scale1	0.04	0.06	495	0.71	0.481
abs:p7_ctry4:scale1	0.01	0.06	495	0.11	0.914

In parallel to Hypothesis #1, the relationship between absorption and secular experience was also significantly positive (see bolded row table above).

For every increase of 1 standard deviation in absorption, participants (on average) increased by 0.36 standard deviations in their secular experiences. In concrete terms, this means that for every increase of 6.51 points on the Absorption scale (range: 0-34), participants (on average) increased by 1.16 points on the Uncommon Events scale (range: 0-18) or 0.73 points on the Extrasensory scale (range: 0-7).

This was actually a (numerically) stronger effect than what we documented for spiritual experiences (b = 0.36 vs. b = 0.21). We did not provide a prediction about whether or not this would be the case in the preregistration.

Parallel to Hypothesis #2

Here are our results, using the Porosity scale (Dulin):

parameter	b	standard error	df	t	p	significant
(Intercept)	0.10	0.04	495	2.36	0.019	*
por	0.45	0.04	495	10.66	<0.001	*
p7_ctry1	0.22	0.07	495	3.21	0.001	*
p7_ctry2	0.04	0.08	495	0.53	0.596
p7_ctry3	-0.21	0.07	495	-3.01	0.003	*
p7_ctry4	0.38	0.09	495	4.06	<0.001	*
scale1	0.10	0.04	495	2.79	0.005	*
por:p7_ctry1	0.24	0.07	495	3.38	<0.001	*
por:p7_ctry2	-0.28	0.07	495	-3.77	<0.001	*
por:p7_ctry3	-0.08	0.08	495	-0.96	0.339
por:p7_ctry4	0.26	0.09	495	2.85	0.005	*
por:scale1	0.14	0.03	495	3.98	<0.001	*
p7_ctry1:scale1	-0.33	0.06	495	-5.79	<0.001	*
p7_ctry2:scale1	-0.20	0.06	495	-3.06	0.002	*
p7_ctry3:scale1	0.15	0.06	495	2.62	0.009	*
p7_ctry4:scale1	0.35	0.08	495	4.59	<0.001	*
por:p7_ctry1:scale1	-0.09	0.06	495	-1.53	0.127
por:p7_ctry2:scale1	-0.12	0.06	495	-1.93	0.054
por:p7_ctry3:scale1	0.08	0.07	495	1.16	0.248
por:p7_ctry4:scale1	0.33	0.07	495	4.45	<0.001	*

And here are our results, using the Martha Story scale:

parameter	b	standard error	df	t	p	significant
(Intercept)	0.09	0.04	495	2.53	0.012	*
mm	0.31	0.04	495	7.92	<0.001	*
p7_ctry1	0.27	0.07	495	3.85	<0.001	*
p7_ctry2	0.05	0.07	495	0.65	0.515
p7_ctry3	-0.30	0.06	495	-4.64	<0.001	*
p7_ctry4	0.10	0.08	495	1.29	0.198
scale1	0.05	0.03	495	1.56	0.118
mm:p7_ctry1	0.29	0.06	495	4.48	<0.001	*
mm:p7_ctry2	-0.12	0.07	495	-1.83	0.067
mm:p7_ctry3	-0.08	0.07	495	-1.07	0.284
mm:p7_ctry4	0.22	0.10	495	2.27	0.023	*
mm:scale1	0.06	0.03	495	1.85	0.065
p7_ctry1:scale1	-0.30	0.06	495	-5.29	<0.001	*
p7_ctry2:scale1	-0.08	0.06	495	-1.39	0.165
p7_ctry3:scale1	0.15	0.05	495	2.80	0.005	*
p7_ctry4:scale1	0.19	0.07	495	2.93	0.004	*
mm:p7_ctry1:scale1	-0.06	0.05	495	-1.09	0.276
mm:p7_ctry2:scale1	-0.15	0.05	495	-2.69	0.007	*
mm:p7_ctry3:scale1	0.05	0.06	495	0.92	0.356
mm:p7_ctry4:scale1	0.24	0.08	495	3.11	0.002	*

In parallel to Hypothesis #2, in both cases, the relationship between porosity and secular experience was significantly positive (see bolded rows in tables above).

For every increase of 1 standard deviation in porosity as indexed by the Porosity (Dulin) scale, participants (on average) increased by 0.45 standard deviations in their secular experiences. In concrete terms, this means that for every increase of 7.47 points on the Porosity scale (range: 0-34), participants (on average) increased by 1.45 points on the Uncommon Events scale (range: 0-18) or 0.91 points on the Extrasensory scale (range: 0-7).

For every increase of 1 standard deviation in porosity as indexed by the Martha Story scale, participants (on average) increased by 0.31 standard deviations in their secular experiences. In concrete terms, this means that for every increase of 10.55 points on the Martha Story scale (range: 0-54), participants (on average) increased by 0.99 points on the Uncommon Events scale (range: 0-18) or 0.62 points on the Extrasensory scale (range: 0-7).

These were actually (numerically) stronger effects than what we documented for spiritual experiences (Porosity: b = 0.45 vs. b = 0.38; Martha Story: b = 0.31 vs. b = 0.26). We did not provide a prediction about whether or not this would be the case in the preregistration.

Parallel to Hypothesis #3

Here are our results:

parameter	b	standard error	df	t	p	significant
US
(Intercept)	0.00	0.07	100	0.00	1
abs	0.46	0.07	100	6.50	<0.001	*
scale1	0.00	0.05	100	0.00	1
abs:scale1	0.00	0.05	100	0.09	0.93
Ghana
(Intercept)	0.00	0.07	98	0.00	1
abs	0.28	0.07	98	4.02	<0.001	*
scale1	0.00	0.07	98	0.00	1
abs:scale1	-0.04	0.07	98	-0.55	0.583
Thailand
(Intercept)	0.00	0.07	109	0.00	1
abs	0.41	0.07	109	6.03	<0.001	*
scale1	0.00	0.06	109	0.00	1
abs:scale1	0.02	0.06	109	0.40	0.693
China
(Intercept)	0.00	0.07	98	0.00	1
abs	0.36	0.07	98	5.12	<0.001	*
scale1	0.00	0.06	98	0.00	1
abs:scale1	-0.07	0.06	98	-1.13	0.263
Vanuatu
(Intercept)	0.00	0.08	90	0.00	1
abs	0.27	0.08	90	3.37	0.001	*
scale1	0.00	0.06	90	0.00	1
abs:scale1	-0.02	0.06	90	-0.37	0.709

In parallel to Hypothesis #3, the relationship between absorption and secular experience was significantly positive in all sites, including Vanuatu (where the relationship with spiritual experience was not significantly positive; see bolded rows in table above). We did not provide a prediction about whether or not this would be the case in the preregistration.

(See other relevant plots under “Parallel to Hypothesis #1,” above.)

Parallel to Hypothesis #4

parameter	b	standard error	df	t	p	significant
US
(Intercept)	0.00	0.07	100	0.00	1
por	0.53	0.07	100	8.03	<0.001	*
scale1	0.00	0.05	100	0.00	1
por:scale1	0.06	0.05	100	1.07	0.286
Ghana
(Intercept)	0.00	0.07	98	0.00	1
por	0.15	0.07	98	2.10	0.038	*
scale1	0.00	0.07	98	0.00	1
por:scale1	0.03	0.07	98	0.41	0.685
Thailand
(Intercept)	0.00	0.07	109	0.00	1
por	0.32	0.07	109	4.37	<0.001	*
scale1	0.00	0.05	109	0.00	1
por:scale1	0.18	0.05	109	3.37	0.001	*
China
(Intercept)	0.00	0.06	98	0.00	1
por	0.46	0.06	98	7.13	<0.001	*
scale1	0.00	0.06	98	0.00	1
por:scale1	0.25	0.06	98	4.31	<0.001	*
Vanuatu
(Intercept)	0.00	0.08	90	0.00	1
por	0.23	0.08	90	2.80	0.006	*
scale1	0.00	0.06	90	0.00	1
por:scale1	-0.05	0.06	90	-0.87	0.384

In parallel to Hypothesis #4, using the Porosity (Dulin) scale, the relationship between porosity and secular experience was significantly positive in all sites, including Vanuatu (where the relationship with spiritual experience was not significantly positive; see bolded rows in table above). We did not provide a prediction about whether or not this would be the case in the preregistration.

And here are our results, using the Martha Story scale:

parameter	b	standard error	df	t	p	significant
US
(Intercept)	0.00	0.07	100	0.00	1
mm	0.52	0.07	100	7.82	<0.001	*
scale1	0.00	0.05	100	0.00	1
mm:scale1	0.02	0.05	100	0.35	0.726
Ghana
(Intercept)	0.00	0.07	98	0.00	1
mm	0.19	0.07	98	2.66	0.009	*
scale1	0.00	0.07	98	0.00	1
mm:scale1	-0.08	0.07	98	-1.18	0.242
Thailand
(Intercept)	0.00	0.07	109	0.00	1
mm	0.24	0.07	109	3.16	0.002	*
scale1	0.00	0.05	109	0.00	1
mm:scale1	0.11	0.05	109	2.02	0.046	*
China
(Intercept)	0.00	0.07	98	0.00	1
mm	0.33	0.07	98	4.58	<0.001	*
scale1	0.00	0.06	98	0.00	1
mm:scale1	0.14	0.06	98	2.30	0.024	*
Vanuatu
(Intercept)	0.00	0.08	90	0.00	1
mm	0.01	0.08	90	0.06	0.95
scale1	0.00	0.06	90	0.00	1
mm:scale1	-0.03	0.06	90	-0.47	0.636

In parallel to Hypothesis #4, using the Martha Story scale, the relationship between porosity and secular experience was significantly positive in all sites except Vanuatu (where the relationship with spiritual experience was also not significantly positive; see bolded rows in table above). We did not provide a prediction about whether or not this would be the case in the preregistration.

(See other relevant plots under “Parallel to Hypothesis #2,” above.)

“Control” scales (Hypothesis #5)

From preregistration: “Second, in order to examine whether the hypothesized relationships between absorption and experience and between porosity and experience can be explained by a lower-level processes—namely, individual or cross-cultural differences in an overall tendency to agree with the questions posed by these scales (i.e., a “yes” bias), we will also examine the relationships between the two “control” scales (the Sense of Control scale, Lachman & Weaver, 1998; and the Need for Cognition scale, Cacioppo, Petty, & Kao, 1993) and our measures of absorption, porosity, and experience. This provides something of a negative control for our primary hypotheses. Our prediction is that any relationships with these control measures will be weaker than the relationships among absorption, porosity, and experience. To assess this, we will standardize scores on each measure, collapsing across field sites, and then conduct a mixed-effects linear regression, using scores on absorption, porosity, and control scales to predict reports of spiritual experiences. Using the “stats” package for R, we will specify this model as follows: experience ~ score * predictor * site + (1 | subject), where “site” is a categorical variable comparing the five field sites (effect-coded), and—critically—“predictor” is a categorical variable with the following orthogonal contrasts: (1) Scales of interest (absorption and the two porosity measures) vs. control scales; (2) absorption vs. porosity measures; (3) porosity as indexed by the Dulin scale vs. porosity as indexed by the Mental Causation scale; (4) Need for Cognition vs. Need for Control (Mastery Subscale). Because of the unusual and complex formulation of this model—considering each of these 5 “predictors” as repeated observations of the variable “score”—we will conduct separate models for the two spiritual experience scales (Daily Spiritual Experiences and Spiritual Events), rather than incorporating both spiritual experiences measures into the same model (as in Hypotheses #1-4). As in Hypotheses #1-4, if we have trouble fitting these models or there are other signs of model fit problems, we will consider the following alternative model: experience ~ score * predictor + (1 | site/subject). We consider these to be confirmatory tests, and in each case our prediction is that the two-way interaction between “score” and the first orthogonal contrast on “predictor” (comparing the predictive power of our scales of interest vs. our two control scales) will be significantly positive (b > 0, t > 2, p < 0.05).”

Although we did not number this hypothesis in the preregistration, I’m going to go ahead and call it Hypothesis #5.

NOTE: Neither of the models we specified in the preregistration were able to run. I think this is due to an error in my thinking (about which of our measures were “repeated measures”). The models reported here are very close to what we preregistered.

Here are our results using the Daily Spiritual Experiences scale:

parameter	b	standard error	t	p	significant
(Intercept)	0.03	0.02	1.95	0.051
score	0.17	0.02	9.61	<0.001	*
predictorint_ctl	0.01	0.01	1.97	0.048	*
predictorabs_pmm	-0.01	0.01	-1.02	0.309
predictorpor_mm	0.00	0.03	0.12	0.906
predictorhthk_wob	0.03	0.03	0.95	0.343
p7_ctry1	-0.37	0.03	-12.07	<0.001	*
p7_ctry2	0.79	0.03	25.22	<0.001	*
p7_ctry3	-0.11	0.04	-3.21	0.001	*
p7_ctry4	-0.96	0.03	-28.68	<0.001	*
score:predictorint_ctl	0.04	0.01	5.57	<0.001	*
score:predictorabs_pmm	-0.04	0.01	-2.76	0.006	*
score:predictorpor_mm	0.08	0.03	2.91	0.004	*
score:predictorhthk_wob	0.02	0.03	0.72	0.473
score:p7_ctry1	0.10	0.03	3.64	<0.001	*
score:p7_ctry2	-0.06	0.03	-1.78	0.075
score:p7_ctry3	0.02	0.04	0.60	0.546
score:p7_ctry4	0.02	0.04	0.53	0.597
predictorint_ctl:p7_ctry1	0.03	0.01	2.12	0.034	*
predictorabs_pmm:p7_ctry1	-0.03	0.03	-1.24	0.214
predictorpor_mm:p7_ctry1	0.01	0.05	0.18	0.855
predictorhthk_wob:p7_ctry1	0.02	0.05	0.38	0.7
predictorint_ctl:p7_ctry2	-0.01	0.01	-1.04	0.298
predictorabs_pmm:p7_ctry2	0.04	0.03	1.35	0.176
predictorpor_mm:p7_ctry2	-0.04	0.05	-0.71	0.478
predictorhthk_wob:p7_ctry2	0.01	0.05	0.30	0.763
predictorint_ctl:p7_ctry3	0.00	0.02	0.22	0.826
predictorabs_pmm:p7_ctry3	0.03	0.03	1.27	0.206
predictorpor_mm:p7_ctry3	0.05	0.05	0.98	0.327
predictorhthk_wob:p7_ctry3	0.01	0.07	0.14	0.885
predictorint_ctl:p7_ctry4	0.02	0.01	1.43	0.153
predictorabs_pmm:p7_ctry4	-0.10	0.03	-3.45	<0.001	*
predictorpor_mm:p7_ctry4	0.10	0.06	1.70	0.089
predictorhthk_wob:p7_ctry4	-0.04	0.05	-0.91	0.361
score:predictorint_ctl:p7_ctry1	0.00	0.01	0.05	0.959
score:predictorabs_pmm:p7_ctry1	-0.05	0.02	-2.00	0.046	*
score:predictorpor_mm:p7_ctry1	0.02	0.05	0.34	0.733
score:predictorhthk_wob:p7_ctry1	-0.07	0.05	-1.46	0.145
score:predictorint_ctl:p7_ctry2	-0.04	0.01	-3.10	0.002	*
score:predictorabs_pmm:p7_ctry2	0.05	0.03	1.75	0.081
score:predictorpor_mm:p7_ctry2	-0.03	0.05	-0.56	0.574
score:predictorhthk_wob:p7_ctry2	-0.08	0.05	-1.64	0.101
score:predictorint_ctl:p7_ctry3	0.01	0.02	0.56	0.574
score:predictorabs_pmm:p7_ctry3	0.04	0.03	1.35	0.177
score:predictorpor_mm:p7_ctry3	0.00	0.05	-0.05	0.958
score:predictorhthk_wob:p7_ctry3	0.04	0.06	0.59	0.558
score:predictorint_ctl:p7_ctry4	0.03	0.02	2.13	0.033	*
score:predictorabs_pmm:p7_ctry4	-0.04	0.03	-1.38	0.166
score:predictorpor_mm:p7_ctry4	-0.01	0.06	-0.08	0.933
score:predictorhthk_wob:p7_ctry4	0.05	0.06	0.81	0.419

And here are our results using the Spiritual Events scale:

parameter	b	standard error	t	p	significant
(Intercept)	0.05	0.02	2.55	0.011	*
score	0.19	0.02	8.68	<0.001	*
predictorint_ctl	0.00	0.01	-0.03	0.978
predictorabs_pmm	-0.01	0.02	-0.39	0.696
predictorpor_mm	0.00	0.03	0.05	0.963
predictorhthk_wob	-0.02	0.03	-0.57	0.57
p7_ctry1	-0.20	0.04	-5.62	<0.001	*
p7_ctry2	0.67	0.04	18.12	<0.001	*
p7_ctry3	-0.26	0.04	-6.17	<0.001	*
p7_ctry4	-0.64	0.04	-16.21	<0.001	*
score:predictorint_ctl	0.06	0.01	6.62	<0.001	*
score:predictorabs_pmm	-0.03	0.02	-2.02	0.043	*
score:predictorpor_mm	0.03	0.03	1.03	0.301
score:predictorhthk_wob	0.08	0.04	2.29	0.022	*
score:p7_ctry1	0.12	0.03	3.42	<0.001	*
score:p7_ctry2	-0.02	0.04	-0.60	0.548
score:p7_ctry3	0.07	0.04	1.72	0.085
score:p7_ctry4	-0.06	0.04	-1.37	0.17
predictorint_ctl:p7_ctry1	0.03	0.02	2.04	0.041	*
predictorabs_pmm:p7_ctry1	-0.05	0.03	-1.77	0.077
predictorpor_mm:p7_ctry1	-0.02	0.06	-0.43	0.664
predictorhthk_wob:p7_ctry1	0.00	0.06	-0.06	0.952
predictorint_ctl:p7_ctry2	-0.04	0.02	-2.55	0.011	*
predictorabs_pmm:p7_ctry2	0.08	0.03	2.51	0.012	*
predictorpor_mm:p7_ctry2	-0.04	0.06	-0.67	0.505
predictorhthk_wob:p7_ctry2	-0.03	0.06	-0.47	0.641
predictorint_ctl:p7_ctry3	0.00	0.02	0.07	0.942
predictorabs_pmm:p7_ctry3	0.03	0.03	1.05	0.296
predictorpor_mm:p7_ctry3	0.04	0.06	0.79	0.428
predictorhthk_wob:p7_ctry3	-0.02	0.08	-0.27	0.786
predictorint_ctl:p7_ctry4	0.02	0.02	1.21	0.226
predictorabs_pmm:p7_ctry4	-0.08	0.04	-2.27	0.023	*
predictorpor_mm:p7_ctry4	0.06	0.07	0.80	0.426
predictorhthk_wob:p7_ctry4	-0.01	0.06	-0.20	0.841
score:predictorint_ctl:p7_ctry1	0.03	0.01	2.44	0.015	*
score:predictorabs_pmm:p7_ctry1	-0.06	0.03	-2.14	0.033	*
score:predictorpor_mm:p7_ctry1	-0.02	0.06	-0.31	0.755
score:predictorhthk_wob:p7_ctry1	-0.04	0.06	-0.78	0.434
score:predictorint_ctl:p7_ctry2	0.03	0.02	1.76	0.078
score:predictorabs_pmm:p7_ctry2	-0.01	0.03	-0.25	0.801
score:predictorpor_mm:p7_ctry2	0.02	0.06	0.38	0.707
score:predictorhthk_wob:p7_ctry2	0.00	0.06	0.06	0.954
score:predictorint_ctl:p7_ctry3	-0.02	0.02	-1.02	0.308
score:predictorabs_pmm:p7_ctry3	0.04	0.03	1.18	0.237
score:predictorpor_mm:p7_ctry3	0.01	0.06	0.24	0.813
score:predictorhthk_wob:p7_ctry3	-0.08	0.07	-1.09	0.276
score:predictorint_ctl:p7_ctry4	0.01	0.02	0.36	0.722
score:predictorabs_pmm:p7_ctry4	-0.01	0.04	-0.38	0.706
score:predictorpor_mm:p7_ctry4	-0.01	0.08	-0.08	0.938
score:predictorhthk_wob:p7_ctry4	0.02	0.07	0.25	0.805

In both cases, we found evidence in favor of Hypothesis #5: As predicted, the interaction between “score” and the first contrast on the “predictor” parameter was significantly positive (see bolded row in tables above), whether indexed spiritual experiences using the DSE or SE scales. In other words, the “predictors of interest” (absorption and the two porosity measures) were better predictors of spiritual experience than our two “control” predictors.

“Control” scales, without reverse coding

The two control scales are the only scales that have reverse-coded items! This is really too bad, because it means that they are not actually good ways to assess a yes-bias.

In an attempt to remedy this, here I re-scoring these scales after dropping reverse-coded items.

Collapsing across sites

Here are our results using the Daily Spiritual Experiences scale:

parameter	b	standard error	t	p	significant
(Intercept)	0.03	0.02	1.94	0.052
score	0.18	0.02	10.57	<0.001	*
predictorint_ctl	0.01	0.01	2.16	0.031	*
predictorabs_pmm	-0.01	0.01	-1.02	0.308
predictorpor_mm	0.00	0.03	0.12	0.906
predictorhthk_wob	0.02	0.02	0.90	0.37
p7_ctry1	-0.35	0.03	-11.84	<0.001	*
p7_ctry2	0.79	0.03	26.04	<0.001	*
p7_ctry3	-0.10	0.03	-3.19	0.001	*
p7_ctry4	-0.96	0.03	-28.80	<0.001	*
score:predictorint_ctl	0.04	0.01	6.03	<0.001	*
score:predictorabs_pmm	-0.04	0.01	-2.76	0.006	*
score:predictorpor_mm	0.08	0.03	2.92	0.004	*
score:predictorhthk_wob	0.08	0.03	3.14	0.002	*
score:p7_ctry1	0.07	0.03	2.57	0.01	*
score:p7_ctry2	-0.06	0.03	-1.95	0.052
score:p7_ctry3	0.04	0.03	1.33	0.183
score:p7_ctry4	0.01	0.04	0.36	0.721
predictorint_ctl:p7_ctry1	0.02	0.01	1.50	0.135
predictorabs_pmm:p7_ctry1	-0.03	0.03	-1.25	0.213
predictorpor_mm:p7_ctry1	0.01	0.05	0.18	0.855
predictorhthk_wob:p7_ctry1	-0.01	0.05	-0.26	0.795
predictorint_ctl:p7_ctry2	-0.02	0.01	-1.24	0.217
predictorabs_pmm:p7_ctry2	0.04	0.03	1.36	0.175
predictorpor_mm:p7_ctry2	-0.04	0.05	-0.71	0.477
predictorhthk_wob:p7_ctry2	-0.03	0.05	-0.59	0.557
predictorint_ctl:p7_ctry3	-0.01	0.01	-0.44	0.663
predictorabs_pmm:p7_ctry3	0.03	0.03	1.27	0.204
predictorpor_mm:p7_ctry3	0.05	0.05	0.98	0.325
predictorhthk_wob:p7_ctry3	0.03	0.05	0.64	0.525
predictorint_ctl:p7_ctry4	0.02	0.01	1.19	0.236
predictorabs_pmm:p7_ctry4	-0.10	0.03	-3.46	<0.001	*
predictorpor_mm:p7_ctry4	0.10	0.06	1.71	0.088
predictorhthk_wob:p7_ctry4	-0.04	0.05	-0.88	0.378
score:predictorint_ctl:p7_ctry1	0.02	0.01	1.42	0.155
score:predictorabs_pmm:p7_ctry1	-0.05	0.02	-2.00	0.045	*
score:predictorpor_mm:p7_ctry1	0.02	0.05	0.34	0.732
score:predictorhthk_wob:p7_ctry1	-0.12	0.05	-2.55	0.011	*
score:predictorint_ctl:p7_ctry2	-0.04	0.01	-3.37	<0.001	*
score:predictorabs_pmm:p7_ctry2	0.05	0.03	1.75	0.08
score:predictorpor_mm:p7_ctry2	-0.03	0.05	-0.56	0.573
score:predictorhthk_wob:p7_ctry2	-0.05	0.04	-1.18	0.24
score:predictorint_ctl:p7_ctry3	0.00	0.01	-0.19	0.847
score:predictorabs_pmm:p7_ctry3	0.04	0.03	1.35	0.176
score:predictorpor_mm:p7_ctry3	0.00	0.05	-0.05	0.957
score:predictorhthk_wob:p7_ctry3	0.00	0.05	-0.01	0.992
score:predictorint_ctl:p7_ctry4	0.04	0.01	2.41	0.016	*
score:predictorabs_pmm:p7_ctry4	-0.04	0.03	-1.39	0.165
score:predictorpor_mm:p7_ctry4	-0.01	0.06	-0.08	0.933
score:predictorhthk_wob:p7_ctry4	0.03	0.06	0.51	0.609

And here are our results using the Spiritual Events scale:

parameter	b	standard error	t	p	significant
(Intercept)	0.03	0.02	1.54	0.123
score	0.20	0.02	10.10	<0.001	*
predictorint_ctl	0.01	0.01	1.35	0.176
predictorabs_pmm	-0.01	0.02	-0.39	0.693
predictorpor_mm	0.00	0.03	0.05	0.962
predictorhthk_wob	-0.01	0.03	-0.29	0.771
p7_ctry1	-0.18	0.04	-5.22	<0.001	*
p7_ctry2	0.66	0.04	18.57	<0.001	*
p7_ctry3	-0.25	0.04	-7.22	<0.001	*
p7_ctry4	-0.62	0.04	-15.94	<0.001	*
score:predictorint_ctl	0.05	0.01	6.67	<0.001	*
score:predictorabs_pmm	-0.03	0.02	-2.04	0.041	*
score:predictorpor_mm	0.03	0.03	1.04	0.297
score:predictorhthk_wob	0.13	0.03	4.46	<0.001	*
score:p7_ctry1	0.10	0.03	3.07	0.002	*
score:p7_ctry2	0.00	0.03	-0.02	0.984
score:p7_ctry3	0.07	0.04	1.71	0.088
score:p7_ctry4	-0.07	0.04	-1.61	0.108
predictorint_ctl:p7_ctry1	0.02	0.01	1.40	0.161
predictorabs_pmm:p7_ctry1	-0.05	0.03	-1.79	0.074
predictorpor_mm:p7_ctry1	-0.02	0.06	-0.44	0.661
predictorhthk_wob:p7_ctry1	-0.04	0.06	-0.73	0.467
predictorint_ctl:p7_ctry2	-0.04	0.01	-2.45	0.014	*
predictorabs_pmm:p7_ctry2	0.08	0.03	2.54	0.011	*
predictorpor_mm:p7_ctry2	-0.04	0.06	-0.67	0.501
predictorhthk_wob:p7_ctry2	-0.04	0.05	-0.69	0.492
predictorint_ctl:p7_ctry3	0.00	0.01	-0.02	0.984
predictorabs_pmm:p7_ctry3	0.03	0.03	1.06	0.291
predictorpor_mm:p7_ctry3	0.04	0.06	0.80	0.423
predictorhthk_wob:p7_ctry3	0.05	0.06	0.85	0.396
predictorint_ctl:p7_ctry4	0.01	0.02	0.53	0.595
predictorabs_pmm:p7_ctry4	-0.08	0.04	-2.29	0.022	*
predictorpor_mm:p7_ctry4	0.06	0.07	0.80	0.422
predictorhthk_wob:p7_ctry4	-0.01	0.05	-0.26	0.796
score:predictorint_ctl:p7_ctry1	0.04	0.01	3.03	0.002	*
score:predictorabs_pmm:p7_ctry1	-0.06	0.03	-2.15	0.031	*
score:predictorpor_mm:p7_ctry1	-0.02	0.05	-0.32	0.753
score:predictorhthk_wob:p7_ctry1	-0.03	0.05	-0.62	0.534
score:predictorint_ctl:p7_ctry2	0.02	0.01	1.17	0.241
score:predictorabs_pmm:p7_ctry2	-0.01	0.03	-0.25	0.8
score:predictorpor_mm:p7_ctry2	0.02	0.06	0.38	0.704
score:predictorhthk_wob:p7_ctry2	0.06	0.05	1.20	0.232
score:predictorint_ctl:p7_ctry3	-0.01	0.02	-0.94	0.348
score:predictorabs_pmm:p7_ctry3	0.04	0.03	1.19	0.233
score:predictorpor_mm:p7_ctry3	0.01	0.06	0.24	0.812
score:predictorhthk_wob:p7_ctry3	-0.11	0.06	-1.80	0.072
score:predictorint_ctl:p7_ctry4	0.01	0.02	0.62	0.534
score:predictorabs_pmm:p7_ctry4	-0.01	0.04	-0.38	0.704
score:predictorpor_mm:p7_ctry4	-0.01	0.08	-0.08	0.937
score:predictorhthk_wob:p7_ctry4	-0.04	0.07	-0.64	0.523

Again, in both cases, we found evidence in favor of Hypothesis #5: As predicted, the interaction between “score” and the first contrast on the “predictor” parameter was significantly positive (see bolded row in tables above), whether indexed spiritual experiences using the DSE or SE scales. In other words, the “predictors of interest” (absorption and the two porosity measures) were better predictors of spiritual experience than our two “control” predictors, even when we only included the non-reverse-coded items from the two “control” predictors.

Removed 15 rows containing non-finite values (stat_smooth).Removed 15 rows containing missing values (geom_point).Removed 15 rows containing non-finite values (stat_smooth).Removed 15 rows containing missing values (geom_point).

Removed 15 rows containing non-finite values (stat_smooth).Removed 15 rows containing missing values (geom_point).Removed 15 rows containing non-finite values (stat_smooth).Removed 15 rows containing missing values (geom_point).

In each site considered alone

Here I conduct a set of parallel analyses, considering each site alone.

Here are our results for the Daily Spiritual Experiences scale:

parameter	b	standard error	t	p	significant
US
(Intercept)	-0.01	0.04	-0.16	0.871
score	0.26	0.04	6.06	<0.001	*
predictorint_ctl	0.00	0.02	0.20	0.842
predictorabs_pmm	0.00	0.04	0.00	1
predictorpor_mm	0.00	0.07	0.00	1
predictorhthk_wob	-0.01	0.07	-0.20	0.844
score:predictorint_ctl	0.06	0.02	3.44	<0.001	*
score:predictorabs_pmm	-0.07	0.04	-1.91	0.057
score:predictorpor_mm	0.07	0.07	1.04	0.298
score:predictorhthk_wob	-0.04	0.07	-0.53	0.598
Ghana
(Intercept)	-0.01	0.04	-0.21	0.836
score	0.19	0.04	4.18	<0.001	*
predictorint_ctl	0.00	0.02	0.25	0.801
predictorabs_pmm	0.00	0.04	0.00	1
predictorpor_mm	0.00	0.07	0.00	1
predictorhthk_wob	0.00	0.07	0.02	0.984
score:predictorint_ctl	-0.01	0.02	-0.39	0.695
score:predictorabs_pmm	0.01	0.04	0.28	0.779
score:predictorpor_mm	0.07	0.07	0.97	0.331
score:predictorhthk_wob	0.07	0.07	0.98	0.326
Thailand
(Intercept)	0.00	0.04	0.00	1
score	0.27	0.04	6.53	<0.001	*
predictorint_ctl	0.00	0.02	0.00	1
predictorabs_pmm	0.00	0.04	0.00	1
predictorpor_mm	0.00	0.06	0.00	1
predictorhthk_wob	0.00	0.06	0.00	1
score:predictorint_ctl	0.05	0.02	2.78	0.006	*
score:predictorabs_pmm	0.01	0.04	0.17	0.863
score:predictorpor_mm	0.06	0.06	0.98	0.327
score:predictorhthk_wob	0.10	0.06	1.51	0.132
China
(Intercept)	0.00	0.04	0.00	1
score	0.23	0.04	5.45	<0.001	*
predictorint_ctl	0.00	0.02	0.00	1
predictorabs_pmm	0.00	0.04	0.00	1
predictorpor_mm	0.00	0.07	0.00	1
predictorhthk_wob	0.00	0.07	0.00	1
score:predictorint_ctl	0.09	0.02	5.03	<0.001	*
score:predictorabs_pmm	-0.08	0.04	-2.10	0.037	*
score:predictorpor_mm	0.10	0.07	1.49	0.137
score:predictorhthk_wob	0.14	0.07	2.12	0.034	*
Vanuatu
(Intercept)	0.00	0.05	-0.04	0.967
score	0.11	0.05	2.42	0.016	*
predictorint_ctl	0.00	0.02	0.05	0.96
predictorabs_pmm	0.00	0.04	0.00	1
predictorpor_mm	0.00	0.07	0.00	1
predictorhthk_wob	0.01	0.07	0.14	0.893
score:predictorint_ctl	0.04	0.02	2.10	0.036	*
score:predictorabs_pmm	-0.02	0.04	-0.42	0.674
score:predictorpor_mm	0.09	0.07	1.22	0.223
score:predictorhthk_wob	0.27	0.07	3.67	<0.001	*

Here are our results for the Daily Spiritual Experiences scale:

parameter	b	standard error	t	p	significant
US
(Intercept)	0.00	0.04	-0.10	0.919
score	0.31	0.04	7.58	<0.001	*
predictorint_ctl	0.00	0.02	0.12	0.901
predictorabs_pmm	0.00	0.04	0.00	1
predictorpor_mm	0.00	0.06	0.00	1
predictorhthk_wob	-0.01	0.07	-0.09	0.925
score:predictorint_ctl	0.10	0.02	5.85	<0.001	*
score:predictorabs_pmm	-0.07	0.04	-1.98	0.048	*
score:predictorpor_mm	-0.02	0.06	-0.33	0.742
score:predictorhthk_wob	0.10	0.07	1.54	0.123
Ghana
(Intercept)	-0.01	0.04	-0.23	0.821
score	0.17	0.04	3.92	<0.001	*
predictorint_ctl	0.00	0.02	0.28	0.782
predictorabs_pmm	0.00	0.04	0.00	1
predictorpor_mm	0.00	0.07	0.00	1
predictorhthk_wob	0.00	0.07	-0.04	0.969
score:predictorint_ctl	0.06	0.02	3.09	0.002	*
score:predictorabs_pmm	-0.04	0.04	-1.03	0.304
score:predictorpor_mm	0.02	0.07	0.31	0.756
score:predictorhthk_wob	0.20	0.07	2.84	0.005	*
Thailand
(Intercept)	0.00	0.04	0.00	1
score	0.31	0.04	7.64	<0.001	*
predictorint_ctl	0.00	0.02	0.00	1
predictorabs_pmm	0.00	0.04	0.00	1
predictorpor_mm	0.00	0.06	0.00	1
predictorhthk_wob	0.00	0.06	0.00	1
score:predictorint_ctl	0.05	0.02	2.74	0.006	*
score:predictorabs_pmm	0.01	0.04	0.37	0.708
score:predictorpor_mm	0.02	0.06	0.34	0.738
score:predictorhthk_wob	0.03	0.06	0.47	0.637
China
(Intercept)	0.00	0.04	0.00	1
score	0.20	0.04	4.59	<0.001	*
predictorint_ctl	0.00	0.02	0.00	1
predictorabs_pmm	0.00	0.04	0.00	1
predictorpor_mm	0.00	0.07	0.00	1
predictorhthk_wob	0.00	0.07	0.00	1
score:predictorint_ctl	0.09	0.02	5.34	<0.001	*
score:predictorabs_pmm	-0.05	0.04	-1.33	0.183
score:predictorpor_mm	0.05	0.07	0.77	0.444
score:predictorhthk_wob	0.14	0.07	2.13	0.033	*
Vanuatu
(Intercept)	-0.01	0.05	-0.31	0.76
score	0.09	0.05	1.87	0.062
predictorint_ctl	0.01	0.02	0.37	0.709
predictorabs_pmm	0.00	0.04	0.00	1
predictorpor_mm	0.00	0.07	0.00	1
predictorhthk_wob	0.01	0.07	0.19	0.851
score:predictorint_ctl	0.00	0.02	0.22	0.824
score:predictorabs_pmm	0.02	0.04	0.55	0.585
score:predictorpor_mm	0.01	0.07	0.12	0.902
score:predictorhthk_wob	0.23	0.07	3.07	0.002	*

In both cases, we found evidence in favor of Hypothesis #5 in almost all sites: In most sites, as predicted the interaction between “score” and the first contrast on the “predictor” parameter was significantly positive (see bolded rows in tables above), whether indexed spiritual experiences using the DSE or SE scales. In other words, the “predictors of interest” (absorption and the two porosity measures) were better predictors of spiritual experience than our two “control” predictors, even when we only included the non-reverse-coded items from the two “control” predictors. The only exceptions to this were Ghana when using the Daily Spiritual Experiences scale (first table above), and Vanuatu when using the Spiritual Events scale (second table above).

Beyond the preregistration

Combined effects of absorption and porosity on spiritual experiences (absorption controlling for porosity & vice versa)

Collapsing across sites

Here I examine a big model, in which both absorption and porosity together predict spiritual experiences. This allows us to look at the effect of absorption controlling for porosity, the effect of porosity controlling for absorption, and the possibility of an interaction between absorption and porosity. I look at porosity as indexed by the Porosity (Dulin) scale and as indexed by the Martha Story separately, and use both the Daily Spiritual Experience and Spiritual Events scales to index spiritual experience.

Here are our results using the Porosity (Dulin) scale:

parameter	b	standard error	df	t	p	significant
(Intercept)	0.07	0.04	485	1.81	0.07
abs	0.14	0.04	485	3.39	<0.001	*
por	0.33	0.04	485	8.07	<0.001	*
p7_ctry1	-0.21	0.06	485	-3.26	0.001	*
p7_ctry2	0.60	0.07	485	8.67	<0.001	*
p7_ctry3	-0.07	0.07	485	-1.10	0.27
p7_ctry4	-0.71	0.10	485	-6.81	<0.001	*
scale1	0.01	0.03	485	0.44	0.661
abs:por	0.04	0.04	485	1.06	0.288
abs:p7_ctry1	0.01	0.06	485	0.10	0.919
abs:p7_ctry2	0.03	0.07	485	0.36	0.72
abs:p7_ctry3	0.13	0.07	485	1.90	0.058
abs:p7_ctry4	0.02	0.10	485	0.17	0.862
por:p7_ctry1	0.20	0.06	485	3.21	0.001	*
por:p7_ctry2	-0.05	0.07	485	-0.79	0.427
por:p7_ctry3	0.01	0.08	485	0.14	0.886
por:p7_ctry4	0.01	0.09	485	0.12	0.905
abs:scale1	-0.04	0.03	485	-1.37	0.172
por:scale1	0.01	0.03	485	0.24	0.811
p7_ctry1:scale1	-0.08	0.04	485	-1.93	0.054
p7_ctry2:scale1	0.09	0.05	485	2.04	0.042	*
p7_ctry3:scale1	0.06	0.04	485	1.39	0.164
p7_ctry4:scale1	-0.09	0.07	485	-1.30	0.194
abs:por:p7_ctry1	-0.03	0.06	485	-0.57	0.572
abs:por:p7_ctry2	-0.12	0.07	485	-1.76	0.079
abs:por:p7_ctry3	0.01	0.08	485	0.07	0.941
abs:por:p7_ctry4	0.07	0.09	485	0.72	0.473
abs:por:scale1	-0.03	0.03	485	-1.34	0.182
abs:p7_ctry1:scale1	-0.02	0.04	485	-0.47	0.641
abs:p7_ctry2:scale1	-0.01	0.05	485	-0.21	0.834
abs:p7_ctry3:scale1	-0.02	0.05	485	-0.50	0.619
abs:p7_ctry4:scale1	-0.02	0.07	485	-0.26	0.798
por:p7_ctry1:scale1	-0.02	0.04	485	-0.37	0.712
por:p7_ctry2:scale1	-0.14	0.04	485	-3.09	0.002	*
por:p7_ctry3:scale1	-0.06	0.05	485	-1.21	0.226
por:p7_ctry4:scale1	0.11	0.06	485	1.75	0.08
abs:por:p7_ctry1:scale1	0.03	0.04	485	0.81	0.419
abs:por:p7_ctry2:scale1	0.09	0.05	485	1.85	0.064
abs:por:p7_ctry3:scale1	-0.09	0.05	485	-1.63	0.103
abs:por:p7_ctry4:scale1	-0.01	0.06	485	-0.11	0.911

And here are our results using the Martha Story scale:

parameter	b	standard error	df	t	p	significant
(Intercept)	0.07	0.03	485	2.12	0.034	*
abs	0.18	0.04	485	5.12	<0.001	*
mm	0.22	0.04	485	6.23	<0.001	*
p7_ctry1	-0.19	0.06	485	-3.07	0.002	*
p7_ctry2	0.70	0.07	485	10.74	<0.001	*
p7_ctry3	-0.16	0.06	485	-2.58	0.01	*
p7_ctry4	-0.85	0.08	485	-11.06	<0.001	*
scale1	0.01	0.02	485	0.66	0.509
abs:mm	0.00	0.04	485	0.05	0.961
abs:p7_ctry1	-0.08	0.06	485	-1.36	0.173
abs:p7_ctry2	0.02	0.07	485	0.29	0.775
abs:p7_ctry3	0.09	0.07	485	1.31	0.192
abs:p7_ctry4	0.03	0.09	485	0.39	0.698
mm:p7_ctry1	0.19	0.06	485	3.28	0.001	*
mm:p7_ctry2	-0.07	0.06	485	-1.13	0.26
mm:p7_ctry3	-0.01	0.07	485	-0.19	0.848
mm:p7_ctry4	0.04	0.09	485	0.42	0.678
abs:scale1	-0.03	0.02	485	-1.54	0.125
mm:scale1	-0.05	0.02	485	-2.21	0.028	*
p7_ctry1:scale1	-0.10	0.04	485	-2.51	0.013	*
p7_ctry2:scale1	0.09	0.04	485	2.16	0.031	*
p7_ctry3:scale1	0.07	0.04	485	1.89	0.06
p7_ctry4:scale1	-0.17	0.05	485	-3.58	<0.001	*
abs:mm:p7_ctry1	-0.01	0.05	485	-0.24	0.808
abs:mm:p7_ctry2	-0.14	0.07	485	-2.15	0.032	*
abs:mm:p7_ctry3	0.06	0.07	485	0.90	0.366
abs:mm:p7_ctry4	0.13	0.10	485	1.29	0.196
abs:mm:scale1	-0.01	0.02	485	-0.53	0.599
abs:p7_ctry1:scale1	-0.01	0.04	485	-0.38	0.705
abs:p7_ctry2:scale1	0.01	0.04	485	0.15	0.884
abs:p7_ctry3:scale1	0.01	0.04	485	0.29	0.769
abs:p7_ctry4:scale1	0.02	0.05	485	0.39	0.7
mm:p7_ctry1:scale1	-0.05	0.04	485	-1.34	0.182
mm:p7_ctry2:scale1	-0.08	0.04	485	-2.13	0.033	*
mm:p7_ctry3:scale1	-0.02	0.04	485	-0.42	0.674
mm:p7_ctry4:scale1	0.09	0.05	485	1.73	0.084
abs:mm:p7_ctry1:scale1	-0.03	0.03	485	-1.03	0.304
abs:mm:p7_ctry2:scale1	0.04	0.04	485	0.99	0.321
abs:mm:p7_ctry3:scale1	-0.07	0.04	485	-1.69	0.091
abs:mm:p7_ctry4:scale1	-0.02	0.06	485	-0.29	0.772

In both cases, we see evidence for absorption and porosity being positively related to spiritual experiences even after controlling for each other: The relationship between absorption and spiritual experience was significantly positive, controlling for porosity; and the effect of porosity (as indexed by either the Dulin scale or the Martha Story scale) and spiritual experience was also significantly positive, controlling for absorption. There was no strong evidence of an interactive effect of absorption and porosity in either case (see bolded rows in tables above).

In each site considered alone

Here I conduct a set of parallel analyses, considering each site alone.

Here are our results using the Porosity (Dulin) scale:

parameter	b	standard error	df	t	p	significant
US
(Intercept)	0.00	0.07	98	-0.04	0.97
abs	0.18	0.08	98	2.33	0.022	*
por	0.50	0.08	98	6.60	<0.001	*
scale1	0.00	0.04	98	0.03	0.973
abs:por	0.01	0.07	98	0.16	0.87
abs:scale1	-0.06	0.04	98	-1.58	0.117
por:scale1	0.00	0.04	98	-0.08	0.939
abs:por:scale1	-0.01	0.04	98	-0.15	0.884
Ghana
(Intercept)	0.02	0.08	96	0.20	0.845
abs	0.12	0.08	96	1.54	0.128
por	0.26	0.08	96	3.18	0.002	*
scale1	-0.01	0.06	96	-0.12	0.908
abs:por	-0.06	0.08	96	-0.85	0.399
abs:scale1	0.02	0.06	96	0.27	0.785
por:scale1	-0.06	0.06	96	-1.07	0.289
abs:por:scale1	0.03	0.06	96	0.50	0.618
Thailand
(Intercept)	-0.01	0.07	107	-0.13	0.898
abs	0.32	0.07	107	4.47	<0.001	*
por	0.34	0.07	107	4.83	<0.001	*
scale1	0.02	0.05	107	0.53	0.594
abs:por	0.04	0.07	107	0.59	0.553
abs:scale1	-0.02	0.05	107	-0.33	0.739
por:scale1	0.01	0.05	107	0.16	0.875
abs:por:scale1	-0.11	0.04	107	-2.46	0.015	*
China
(Intercept)	-0.04	0.08	96	-0.50	0.617
abs	0.12	0.08	96	1.59	0.116
por	0.48	0.08	96	6.08	<0.001	*
scale1	0.02	0.05	96	0.40	0.691
abs:por	0.12	0.07	96	1.68	0.097
abs:scale1	-0.05	0.05	96	-1.05	0.297
por:scale1	0.08	0.05	96	1.61	0.111
abs:por:scale1	-0.06	0.04	96	-1.33	0.187
Vanuatu
(Intercept)	-0.02	0.09	88	-0.25	0.8
abs	0.11	0.09	88	1.30	0.197
por	0.17	0.09	88	2.01	0.047	*
scale1	0.01	0.06	88	0.11	0.912
abs:por	0.10	0.08	88	1.17	0.246
abs:scale1	-0.02	0.06	88	-0.27	0.789
por:scale1	0.11	0.06	88	1.78	0.079
abs:por:scale1	-0.03	0.06	88	-0.51	0.613

And here are our results using the Martha Story scale:

parameter	b	standard error	df	t	p	significant
US
(Intercept)	0.01	0.08	98	0.06	0.949
abs	0.13	0.08	98	1.61	0.111
mm	0.45	0.08	98	5.37	<0.001	*
scale1	0.02	0.04	98	0.54	0.589
abs:mm	-0.01	0.08	98	-0.19	0.852
abs:scale1	-0.04	0.04	98	-0.87	0.387
mm:scale1	-0.10	0.04	98	-2.43	0.017	*
abs:mm:scale1	-0.06	0.04	98	-1.57	0.119
Ghana
(Intercept)	0.04	0.08	96	0.46	0.646
abs	0.13	0.08	96	1.61	0.111
mm	0.15	0.08	96	1.88	0.064
scale1	0.00	0.06	96	0.04	0.966
abs:mm	-0.15	0.08	96	-1.80	0.076
abs:scale1	0.02	0.06	96	0.40	0.689
mm:scale1	-0.11	0.06	96	-1.92	0.058
abs:mm:scale1	-0.01	0.06	96	-0.17	0.869
Thailand
(Intercept)	-0.02	0.07	107	-0.26	0.797
abs	0.32	0.08	107	4.26	<0.001	*
mm	0.23	0.08	107	3.05	0.003	*
scale1	0.03	0.05	107	0.56	0.575
abs:mm	0.07	0.07	107	0.98	0.329
abs:scale1	0.00	0.05	107	-0.10	0.918
mm:scale1	-0.03	0.05	107	-0.62	0.54
abs:mm:scale1	-0.09	0.04	107	-2.14	0.035	*
China
(Intercept)	-0.01	0.08	96	-0.17	0.862
abs	0.24	0.08	96	2.97	0.004	*
mm	0.37	0.08	96	4.70	<0.001	*
scale1	0.00	0.04	96	0.10	0.917
abs:mm	0.14	0.08	96	1.70	0.092
abs:scale1	-0.02	0.04	96	-0.50	0.619
mm:scale1	0.00	0.04	96	0.00	1
abs:mm:scale1	-0.05	0.05	96	-1.02	0.311
Vanuatu
(Intercept)	0.00	0.09	88	0.04	0.97
abs	0.14	0.09	88	1.63	0.106
mm	0.07	0.09	88	0.82	0.414
scale1	-0.01	0.06	88	-0.14	0.885
abs:mm	-0.03	0.08	88	-0.32	0.752
abs:scale1	0.00	0.06	88	-0.03	0.974
mm:scale1	0.04	0.06	88	0.60	0.55
abs:mm:scale1	0.07	0.06	88	1.22	0.225

In both cases, we see some evidence for porosity “swamping out” the effect of absorption in some sites, including Ghana, China, and Vanuatu when porosity is indexed by the Porosity (Dulin) scale; and the US, Ghana, and Vanuatu when porosity is indexed by the Martha Story scale (in fact, in this second set of analyses, neither Absorption nor Martha Story ends up being a significant predictor of spiritual experience in Ghana or Vanuatu).

Absorption as primarily individual difference, porosity as primarily cultural difference?

In our theorizing, we often assume that absorption is primarily an individual difference, while porosity is primarily a cultural difference—or at least that absorption is more of an individual difference than is porosity. Do we see any evidence of this in this dataset?

Do individuals vary more in absorption than porosity?

The first way I thought of to approach this question is to look at the degree to which individuals vary in each of these three scales (Absorption, Porosity (Dulin), and Martha Story). We would predict that absorption would vary more across individuals than either of the porosity indices (after resclaing all threee scales to have the same range, within each site).

[1] "Absorption vs. Porosity (Dulin)"
[1] "Overall"

    Bartlett test of homogeneity of variances

data:  score by scale
Bartlett's K-squared = 0.16827, df = 1, p-value = 0.6817

[1] "US"

    Bartlett test of homogeneity of variances

data:  score by scale
Bartlett's K-squared = 0.35391, df = 1, p-value = 0.5519

[1] "Ghana"

    Bartlett test of homogeneity of variances

data:  score by scale
Bartlett's K-squared = 3.2876, df = 1, p-value = 0.0698

[1] "Thailand"

    Bartlett test of homogeneity of variances

data:  score by scale
Bartlett's K-squared = 0.46254, df = 1, p-value = 0.4964

[1] "China"

    Bartlett test of homogeneity of variances

data:  score by scale
Bartlett's K-squared = 0.00084535, df = 1, p-value = 0.9768

[1] "Vanuatu"

    Bartlett test of homogeneity of variances

data:  score by scale
Bartlett's K-squared = 9.996, df = 1, p-value = 0.001569

[1] "Absorption vs. Martha Story"
[1] "Overall"

    Bartlett test of homogeneity of variances

data:  score by scale
Bartlett's K-squared = 0.16827, df = 1, p-value = 0.6817

[1] "US"

    Bartlett test of homogeneity of variances

data:  score by scale
Bartlett's K-squared = 0.63721, df = 1, p-value = 0.4247

[1] "Ghana"

    Bartlett test of homogeneity of variances

data:  score by scale
Bartlett's K-squared = 3.8497, df = 1, p-value = 0.04976

[1] "Thailand"

    Bartlett test of homogeneity of variances

data:  score by scale
Bartlett's K-squared = 1.089, df = 1, p-value = 0.2967

[1] "China"

    Bartlett test of homogeneity of variances

data:  score by scale
Bartlett's K-squared = 3.9359, df = 1, p-value = 0.04726

[1] "Vanuatu"

    Bartlett test of homogeneity of variances

data:  score by scale
Bartlett's K-squared = 6.5318, df = 1, p-value = 0.0106

[1] "Dulin vs. Martha Story"
[1] "Overall"

    Bartlett test of homogeneity of variances

data:  score by scale
Bartlett's K-squared = 0.16827, df = 1, p-value = 0.6817

[1] "US"

    Bartlett test of homogeneity of variances

data:  score by scale
Bartlett's K-squared = 0.041448, df = 1, p-value = 0.8387

[1] "Ghana"

    Bartlett test of homogeneity of variances

data:  score by scale
Bartlett's K-squared = 0.022569, df = 1, p-value = 0.8806

[1] "Thailand"

    Bartlett test of homogeneity of variances

data:  score by scale
Bartlett's K-squared = 0.13252, df = 1, p-value = 0.7158

[1] "China"

    Bartlett test of homogeneity of variances

data:  score by scale
Bartlett's K-squared = 4.051, df = 1, p-value = 0.04415

[1] "Vanuatu"

    Bartlett test of homogeneity of variances

data:  score by scale
Bartlett's K-squared = 0.38374, df = 1, p-value = 0.5356

According to Bartlett tests of homogeneity of variances, the only site where we see the pattern that we’d expect is Vanuatu: Individuals differed significantly more in their Absorption scores than in their Porosity (Dulin) scores or Martha Story scores. In China, individuals differed significantly more in their Absorption scores than in their Martha Story scores, but variance was equal across the Absorption and Porosity scores. In Ghana, we actually saw the reverse pattern to what we’d expect: Individuals differed significantly less in their Absorption scores than in their Martha Story scores (and the differences between Absorption and Porosity (Dulin) also went in the same directly, although it was not significant).

My take-away here is that there is no strong evidence from these analyses of homogeneity of variance that individuals (within a site) varied more in absorption than in porosity.

Do site differences account for more of the variance in porosity than in absorption?

The second way I thought of to approach this question is to look at the amount of variance explained by site differences for in each of these three scales (Absorption, Porosity (Dulin), and Martha Story). We would predict that site differences would explain more of the variance in porosity than in absorption (no rescaling necessary here).

Numerically speaking, this appears to be true: In a set of linear models, site differences explain only 9% of the variance in Absorption scores, vs. 41% of the variance in Porosity (Dulin) scores and 22% of the variance in Martha Story scores.

Here, I’ll try to compare these statistically (using scores for each scale that have been standardized, collapsing across sites):

parameter	b	standard error	df	t	p	significant
(Intercept)	0.02	0.03	500	0.62	0.537
p7_ctrynonUS_US	0.07	0.01	500	4.79	<0.001	*
p7_ctryGHVT_THCH	0.41	0.03	500	12.63	<0.001	*
p7_ctryGH_VT	-0.12	0.05	500	-2.45	0.015	*
p7_ctryTH_CH	0.00	0.05	500	-0.01	0.991
scalepmm_abs	0.00	0.01	1000	0.10	0.919
scalepor_mm	0.00	0.02	1000	0.22	0.823
p7_ctrynonUS_US:scalepmm_abs	0.02	0.01	1000	3.69	<0.001	*
p7_ctryGHVT_THCH:scalepmm_abs	0.15	0.01	1000	10.71	<0.001	*
p7_ctryGH_VT:scalepmm_abs	0.01	0.02	1000	0.72	0.473
p7_ctryTH_CH:scalepmm_abs	-0.21	0.02	1000	-10.38	<0.001	*
p7_ctrynonUS_US:scalepor_mm	-0.01	0.01	1000	-1.03	0.305
p7_ctryGHVT_THCH:scalepor_mm	0.11	0.02	1000	4.35	<0.001	*
p7_ctryGH_VT:scalepor_mm	-0.08	0.04	1000	-2.32	0.021	*
p7_ctryTH_CH:scalepor_mm	-0.01	0.03	1000	-0.43	0.667

This analysis appears to be consistent: I think that the interactions that I have bolded in the table above suggest that the differences across sites (here, coded as differences between (a) non-US and US sites, (b) Ghana and Vanuatu vs. Thailand and China, (c) Ghana vs. Vanuatu, and (d) Thailand vs. China) were generally exaggerated in the two porosity scales compared to the Absorption scale, and these seems to be true from the following plot:

My take-away here is that there is pretty strong evidence from these resression analyses that site differences accounted for more of the variance in porosity than in absorption.

Religiosity

Collapsing across sites

Do the relationships between absorption and porosity, on the one hand, and spiritual experience, on the other hand, still hold even after controlling for participants’ self-reported religiosity?

parameter	b	standard error	df	t	p	significant
(Intercept)	0.05	0.05	461	0.90	0.366
abs	0.10	0.05	461	1.82	0.069
por	0.18	0.06	461	3.15	0.002	*
relig_std	0.33	0.06	461	5.10	<0.001	*
scale1	0.02	0.04	461	0.51	0.608
p7_ctry1	-0.22	0.07	461	-2.98	0.003	*
p7_ctry2	0.46	0.08	461	5.78	<0.001	*
p7_ctry3	-0.16	0.08	461	-2.15	0.032	*
p7_ctry4	-0.51	0.14	461	-3.59	<0.001	*
abs:por	0.02	0.06	461	0.38	0.706
abs:relig_std	0.09	0.07	461	1.35	0.177
por:relig_std	0.16	0.06	461	2.72	0.007	*
abs:scale1	-0.04	0.04	461	-1.03	0.303
por:scale1	0.02	0.04	461	0.56	0.573
relig_std:scale1	0.00	0.05	461	0.05	0.963
abs:p7_ctry1	-0.01	0.07	461	-0.15	0.88
abs:p7_ctry2	0.10	0.09	461	1.21	0.226
abs:p7_ctry3	0.13	0.08	461	1.69	0.093
abs:p7_ctry4	0.09	0.13	461	0.65	0.515
por:p7_ctry1	0.10	0.08	461	1.26	0.209
por:p7_ctry2	-0.06	0.09	461	-0.64	0.521
por:p7_ctry3	-0.03	0.10	461	-0.29	0.775
por:p7_ctry4	0.15	0.16	461	0.95	0.344
relig_std:p7_ctry1	0.22	0.08	461	2.68	0.008	*
relig_std:p7_ctry2	-0.01	0.09	461	-0.15	0.882
relig_std:p7_ctry3	0.11	0.11	461	0.93	0.354
relig_std:p7_ctry4	-0.05	0.14	461	-0.34	0.737
scale1:p7_ctry1	-0.08	0.05	461	-1.44	0.151
scale1:p7_ctry2	0.12	0.06	461	2.12	0.034	*
scale1:p7_ctry3	0.09	0.05	461	1.59	0.113
scale1:p7_ctry4	-0.13	0.10	461	-1.23	0.219
abs:por:relig_std	-0.03	0.06	461	-0.49	0.626
abs:por:scale1	-0.03	0.04	461	-0.79	0.429
abs:relig_std:scale1	-0.01	0.05	461	-0.12	0.903
por:relig_std:scale1	-0.09	0.04	461	-1.96	0.05
abs:por:p7_ctry1	-0.17	0.08	461	-2.25	0.025	*
abs:por:p7_ctry2	-0.17	0.10	461	-1.82	0.07
abs:por:p7_ctry3	-0.07	0.09	461	-0.75	0.455
abs:por:p7_ctry4	0.19	0.15	461	1.24	0.215
abs:relig_std:p7_ctry1	-0.11	0.08	461	-1.33	0.185
abs:relig_std:p7_ctry2	-0.15	0.10	461	-1.44	0.151
abs:relig_std:p7_ctry3	-0.05	0.13	461	-0.41	0.684
abs:relig_std:p7_ctry4	0.03	0.13	461	0.22	0.825
por:relig_std:p7_ctry1	0.03	0.07	461	0.42	0.673
por:relig_std:p7_ctry2	-0.07	0.09	461	-0.78	0.434
por:relig_std:p7_ctry3	0.06	0.13	461	0.46	0.644
por:relig_std:p7_ctry4	-0.05	0.12	461	-0.37	0.712
abs:scale1:p7_ctry1	-0.04	0.05	461	-0.80	0.424
abs:scale1:p7_ctry2	-0.08	0.06	461	-1.23	0.219
abs:scale1:p7_ctry3	0.01	0.05	461	0.22	0.826
abs:scale1:p7_ctry4	0.01	0.10	461	0.11	0.913
por:scale1:p7_ctry1	-0.10	0.06	461	-1.73	0.084
por:scale1:p7_ctry2	-0.08	0.06	461	-1.22	0.224
por:scale1:p7_ctry3	-0.02	0.07	461	-0.34	0.733
por:scale1:p7_ctry4	0.07	0.11	461	0.59	0.556
relig_std:scale1:p7_ctry1	0.01	0.06	461	0.18	0.86
relig_std:scale1:p7_ctry2	-0.02	0.07	461	-0.32	0.747
relig_std:scale1:p7_ctry3	-0.03	0.08	461	-0.33	0.742
relig_std:scale1:p7_ctry4	-0.03	0.10	461	-0.36	0.72
abs:por:relig_std:scale1	0.02	0.04	461	0.35	0.727
abs:por:relig_std:p7_ctry1	-0.05	0.07	461	-0.67	0.504
abs:por:relig_std:p7_ctry2	0.08	0.09	461	0.83	0.405
abs:por:relig_std:p7_ctry3	0.06	0.14	461	0.42	0.672
abs:por:relig_std:p7_ctry4	0.18	0.11	461	1.56	0.12
abs:por:scale1:p7_ctry1	0.12	0.05	461	2.19	0.029	*
abs:por:scale1:p7_ctry2	0.06	0.07	461	0.92	0.359
abs:por:scale1:p7_ctry3	-0.03	0.07	461	-0.44	0.658
abs:por:scale1:p7_ctry4	0.02	0.11	461	0.19	0.848
abs:relig_std:scale1:p7_ctry1	-0.04	0.06	461	-0.67	0.501
abs:relig_std:scale1:p7_ctry2	0.12	0.07	461	1.65	0.099
abs:relig_std:scale1:p7_ctry3	-0.09	0.09	461	-1.04	0.301
abs:relig_std:scale1:p7_ctry4	0.04	0.10	461	0.40	0.691
por:relig_std:scale1:p7_ctry1	0.01	0.05	461	0.20	0.845
por:relig_std:scale1:p7_ctry2	-0.01	0.06	461	-0.23	0.82
por:relig_std:scale1:p7_ctry3	-0.03	0.10	461	-0.29	0.772
por:relig_std:scale1:p7_ctry4	0.06	0.09	461	0.71	0.478
abs:por:relig_std:scale1:p7_ctry1	0.05	0.05	461	1.00	0.318
abs:por:relig_std:scale1:p7_ctry2	-0.02	0.07	461	-0.24	0.809
abs:por:relig_std:scale1:p7_ctry3	-0.23	0.10	461	-2.35	0.019	*
abs:por:relig_std:scale1:p7_ctry4	0.01	0.08	461	0.17	0.861

parameter	b	standard error	df	t	p	significant
(Intercept)	0.04	0.04	461	1.12	0.264
abs	0.19	0.04	461	4.48	<0.001	*
mm	0.14	0.05	461	2.95	0.003	*
relig_std	0.36	0.05	461	7.38	<0.001	*
scale1	0.04	0.03	461	1.49	0.138
p7_ctry1	-0.18	0.06	461	-2.87	0.004	*
p7_ctry2	0.51	0.07	461	7.01	<0.001	*
p7_ctry3	-0.18	0.06	461	-3.00	0.003	*
p7_ctry4	-0.60	0.10	461	-6.29	<0.001	*
abs:mm	-0.09	0.05	461	-1.58	0.114
abs:relig_std	0.02	0.05	461	0.44	0.664
mm:relig_std	0.09	0.05	461	1.93	0.054
abs:scale1	-0.06	0.03	461	-1.95	0.052
mm:scale1	-0.04	0.03	461	-1.19	0.233
relig_std:scale1	0.00	0.03	461	0.06	0.95
abs:p7_ctry1	-0.15	0.06	461	-2.50	0.013	*
abs:p7_ctry2	0.05	0.08	461	0.61	0.544
abs:p7_ctry3	0.06	0.07	461	0.93	0.352
abs:p7_ctry4	0.06	0.10	461	0.62	0.536
mm:p7_ctry1	0.12	0.07	461	1.76	0.08
mm:p7_ctry2	-0.06	0.07	461	-0.83	0.406
mm:p7_ctry3	-0.05	0.08	461	-0.62	0.537
mm:p7_ctry4	0.19	0.14	461	1.36	0.175
relig_std:p7_ctry1	0.12	0.07	461	1.89	0.059
relig_std:p7_ctry2	-0.02	0.08	461	-0.20	0.839
relig_std:p7_ctry3	0.09	0.11	461	0.82	0.41
relig_std:p7_ctry4	-0.07	0.09	461	-0.75	0.457
scale1:p7_ctry1	-0.10	0.04	461	-2.19	0.029	*
scale1:p7_ctry2	0.12	0.05	461	2.27	0.024	*
scale1:p7_ctry3	0.06	0.04	461	1.40	0.162
scale1:p7_ctry4	-0.19	0.07	461	-2.79	0.005	*
abs:mm:relig_std	0.09	0.05	461	1.69	0.092
abs:mm:scale1	-0.04	0.04	461	-1.04	0.298
abs:relig_std:scale1	0.03	0.04	461	0.77	0.44
mm:relig_std:scale1	-0.07	0.03	461	-1.96	0.05
abs:mm:p7_ctry1	0.06	0.07	461	0.94	0.349
abs:mm:p7_ctry2	-0.11	0.08	461	-1.35	0.177
abs:mm:p7_ctry3	0.08	0.08	461	1.02	0.31
abs:mm:p7_ctry4	0.17	0.16	461	1.09	0.276
abs:relig_std:p7_ctry1	-0.02	0.07	461	-0.24	0.811
abs:relig_std:p7_ctry2	-0.10	0.09	461	-1.05	0.295
abs:relig_std:p7_ctry3	0.07	0.13	461	0.58	0.56
abs:relig_std:p7_ctry4	0.10	0.11	461	0.95	0.341
mm:relig_std:p7_ctry1	0.02	0.06	461	0.26	0.794
mm:relig_std:p7_ctry2	-0.03	0.07	461	-0.39	0.693
mm:relig_std:p7_ctry3	-0.03	0.10	461	-0.35	0.723
mm:relig_std:p7_ctry4	0.07	0.10	461	0.69	0.488
abs:scale1:p7_ctry1	0.00	0.04	461	-0.02	0.987
abs:scale1:p7_ctry2	-0.04	0.06	461	-0.68	0.494
abs:scale1:p7_ctry3	0.05	0.05	461	1.06	0.292
abs:scale1:p7_ctry4	0.04	0.07	461	0.57	0.569
mm:scale1:p7_ctry1	-0.10	0.05	461	-2.22	0.027	*
mm:scale1:p7_ctry2	-0.02	0.05	461	-0.47	0.64
mm:scale1:p7_ctry3	0.02	0.06	461	0.34	0.737
mm:scale1:p7_ctry4	0.05	0.10	461	0.53	0.596
relig_std:scale1:p7_ctry1	0.05	0.05	461	1.08	0.28
relig_std:scale1:p7_ctry2	-0.05	0.06	461	-0.81	0.417
relig_std:scale1:p7_ctry3	0.01	0.08	461	0.08	0.934
relig_std:scale1:p7_ctry4	0.01	0.06	461	0.12	0.905
abs:mm:relig_std:scale1	0.02	0.04	461	0.54	0.589
abs:mm:relig_std:p7_ctry1	-0.06	0.06	461	-0.90	0.37
abs:mm:relig_std:p7_ctry2	-0.01	0.08	461	-0.12	0.905
abs:mm:relig_std:p7_ctry3	-0.06	0.11	461	-0.55	0.58
abs:mm:relig_std:p7_ctry4	-0.03	0.12	461	-0.25	0.803
abs:mm:scale1:p7_ctry1	0.02	0.05	461	0.34	0.736
abs:mm:scale1:p7_ctry2	0.03	0.06	461	0.59	0.553
abs:mm:scale1:p7_ctry3	0.00	0.06	461	-0.04	0.967
abs:mm:scale1:p7_ctry4	-0.02	0.11	461	-0.19	0.848
abs:relig_std:scale1:p7_ctry1	-0.03	0.05	461	-0.72	0.472
abs:relig_std:scale1:p7_ctry2	0.07	0.06	461	1.07	0.285
abs:relig_std:scale1:p7_ctry3	-0.05	0.09	461	-0.51	0.612
abs:relig_std:scale1:p7_ctry4	-0.04	0.07	461	-0.51	0.61
mm:relig_std:scale1:p7_ctry1	0.00	0.04	461	0.03	0.979
mm:relig_std:scale1:p7_ctry2	-0.05	0.05	461	-0.95	0.345
mm:relig_std:scale1:p7_ctry3	-0.02	0.07	461	-0.30	0.765
mm:relig_std:scale1:p7_ctry4	0.04	0.07	461	0.54	0.589
abs:mm:relig_std:scale1:p7_ctry1	0.02	0.04	461	0.51	0.608
abs:mm:relig_std:scale1:p7_ctry2	0.03	0.06	461	0.53	0.594
abs:mm:relig_std:scale1:p7_ctry3	-0.11	0.08	461	-1.45	0.148
abs:mm:relig_std:scale1:p7_ctry4	-0.05	0.08	461	-0.64	0.522

Following other analysis, I examined this separately using Porosity (Dulin) and Martha Story as indices of porosity. In both cases, religiosity was a significant positive predictor of spiritual experience (as we’d expect). But, importantly, the main effects of porosity remained significantly positive in both analyses, even after controlling for participants’ self-reported religiosity (see bolded rows in tables above). The main effect of absorption remained significantly positive when porosity was indexed by the Martha Story (second table), but not when porosity was indexed by the Porosity (Dulin) scale (although some would consider this effect to be “marginally significant”).

My take-away here is that religiosity does not fully “explain away” the observed relationships between absorption and porosity, on the one hand, and spiritual experience, on the other.

In each site considered alone

Here I conduct a set of parallel analyses, considering each site alone.

parameter	b	standard error	df	t	p	significant
US
(Intercept)	0.04	0.07	93	0.50	0.617
abs	0.16	0.07	93	2.19	0.031	*
por	0.28	0.08	93	3.36	0.001	*
relig_std	0.49	0.07	93	7.12	<0.001	*
scale1	0.04	0.04	93	0.92	0.361
abs:por	-0.17	0.08	93	-2.01	0.047	*
abs:relig_std	0.01	0.07	93	0.15	0.884
por:relig_std	0.19	0.06	93	3.10	0.003	*
abs:scale1	-0.13	0.05	93	-2.80	0.006	*
por:scale1	-0.08	0.05	93	-1.55	0.125
relig_std:scale1	0.05	0.04	93	1.07	0.288
abs:por:relig_std	-0.09	0.06	93	-1.41	0.161
abs:por:scale1	0.10	0.05	93	1.89	0.061
abs:relig_std:scale1	-0.08	0.04	93	-1.85	0.067
por:relig_std:scale1	-0.07	0.04	93	-1.96	0.053
abs:por:relig_std:scale1	0.08	0.04	93	1.98	0.051
Ghana
(Intercept)	-0.28	0.11	91	-2.51	0.014	*
abs	0.09	0.11	91	0.85	0.398
por	0.12	0.10	91	1.17	0.244
relig_std	0.44	0.12	91	3.73	<0.001	*
scale1	0.01	0.09	91	0.10	0.917
abs:por	-0.14	0.10	91	-1.33	0.188
abs:relig_std	0.01	0.11	91	0.04	0.965
por:relig_std	0.08	0.09	91	0.83	0.408
abs:scale1	-0.07	0.08	91	-0.87	0.386
por:scale1	-0.02	0.08	91	-0.28	0.782
relig_std:scale1	0.03	0.09	91	0.28	0.778
abs:por:relig_std	0.04	0.09	91	0.48	0.633
abs:por:scale1	-0.01	0.08	91	-0.12	0.906
abs:relig_std:scale1	0.12	0.09	91	1.31	0.194
por:relig_std:scale1	-0.08	0.07	91	-1.10	0.273
abs:por:relig_std:scale1	0.01	0.07	91	0.16	0.877
Thailand
(Intercept)	-0.08	0.07	102	-1.12	0.266
abs	0.30	0.07	102	4.56	<0.001	*
por	0.19	0.07	102	2.59	0.011	*
relig_std	0.51	0.13	102	3.98	<0.001	*
scale1	0.03	0.05	102	0.71	0.48
abs:por	-0.04	0.07	102	-0.62	0.539
abs:relig_std	0.04	0.15	102	0.24	0.811
por:relig_std	0.23	0.13	102	1.73	0.087
abs:scale1	0.00	0.05	102	-0.04	0.968
por:scale1	0.04	0.05	102	0.70	0.487
relig_std:scale1	0.05	0.09	102	0.56	0.577
abs:por:relig_std	0.02	0.12	102	0.18	0.861
abs:por:scale1	-0.06	0.05	102	-1.15	0.255
abs:relig_std:scale1	-0.03	0.11	102	-0.27	0.784
por:relig_std:scale1	0.00	0.09	102	-0.05	0.96
abs:por:relig_std:scale1	-0.20	0.08	102	-2.30	0.023	*
China
(Intercept)	0.27	0.12	92	2.19	0.031	*
abs	0.03	0.14	92	0.25	0.804
por	0.52	0.13	92	4.03	<0.001	*
relig_std	0.38	0.10	92	3.80	<0.001	*
scale1	-0.04	0.09	92	-0.42	0.672
abs:por	0.23	0.12	92	1.87	0.065
abs:relig_std	0.00	0.10	92	0.04	0.968
por:relig_std	0.23	0.10	92	2.36	0.02	*
abs:scale1	-0.04	0.10	92	-0.38	0.707
por:scale1	0.05	0.09	92	0.61	0.541
relig_std:scale1	-0.06	0.07	92	-0.94	0.352
abs:por:relig_std	0.15	0.08	92	1.87	0.064
abs:por:scale1	-0.04	0.08	92	-0.45	0.657
abs:relig_std:scale1	0.01	0.07	92	0.18	0.855
por:relig_std:scale1	-0.04	0.07	92	-0.64	0.527
abs:por:relig_std:scale1	0.01	0.06	92	0.25	0.807
Vanuatu
(Intercept)	-0.22	0.12	83	-1.88	0.063
abs	0.01	0.14	83	0.04	0.972
por	0.08	0.12	83	0.66	0.514
relig_std	0.38	0.15	83	2.50	0.014	*
scale1	-0.04	0.08	83	-0.45	0.652
abs:por	0.18	0.14	83	1.30	0.199
abs:relig_std	0.15	0.17	83	0.89	0.377
por:relig_std	0.07	0.15	83	0.49	0.626
abs:scale1	-0.21	0.10	83	-2.03	0.046	*
por:scale1	0.09	0.09	83	1.07	0.287
relig_std:scale1	0.07	0.11	83	0.67	0.503
abs:por:relig_std	-0.22	0.18	83	-1.24	0.218
abs:por:scale1	-0.14	0.10	83	-1.44	0.154
abs:relig_std:scale1	0.25	0.12	83	2.07	0.041	*
por:relig_std:scale1	-0.04	0.11	83	-0.34	0.737
abs:por:relig_std:scale1	0.13	0.13	83	1.00	0.32

And here are our results using the Martha Story scale:

parameter	b	standard error	df	t	p	significant
US
(Intercept)	0.06	0.07	93	0.78	0.44
abs	0.06	0.07	93	0.90	0.371
mm	0.28	0.07	93	3.92	<0.001	*
relig_std	0.44	0.06	93	7.20	<0.001	*
scale1	0.07	0.04	93	1.67	0.098
abs:mm	-0.03	0.07	93	-0.47	0.641
abs:relig_std	-0.01	0.06	93	-0.13	0.898
mm:relig_std	0.11	0.06	93	2.01	0.047	*
abs:scale1	-0.06	0.04	93	-1.43	0.155
mm:scale1	-0.15	0.04	93	-3.63	<0.001	*
relig_std:scale1	0.08	0.04	93	2.39	0.019	*
abs:mm:relig_std	0.05	0.06	93	0.71	0.48
abs:mm:scale1	-0.03	0.04	93	-0.79	0.432
abs:relig_std:scale1	-0.03	0.03	93	-0.77	0.446
mm:relig_std:scale1	-0.07	0.03	93	-2.22	0.029	*
abs:mm:relig_std:scale1	0.06	0.04	93	1.54	0.126
Ghana
(Intercept)	-0.25	0.11	91	-2.40	0.018	*
abs	0.11	0.10	91	1.08	0.285
mm	0.09	0.09	91	0.96	0.339
relig_std	0.43	0.11	91	3.94	<0.001	*
scale1	0.03	0.08	91	0.37	0.713
abs:mm	-0.22	0.10	91	-2.30	0.024	*
abs:relig_std	0.01	0.11	91	0.07	0.945
mm:relig_std	0.04	0.09	91	0.50	0.62
abs:scale1	-0.07	0.08	91	-0.94	0.349
mm:scale1	-0.05	0.07	91	-0.76	0.448
relig_std:scale1	-0.01	0.08	91	-0.11	0.916
abs:mm:relig_std	0.11	0.09	91	1.18	0.242
abs:mm:scale1	-0.07	0.07	91	-0.89	0.376
abs:relig_std:scale1	0.13	0.08	91	1.62	0.109
mm:relig_std:scale1	-0.12	0.07	91	-1.79	0.076
abs:mm:relig_std:scale1	0.08	0.07	91	1.18	0.241
Thailand
(Intercept)	-0.06	0.07	102	-0.87	0.387
abs	0.31	0.07	102	4.37	<0.001	*
mm	0.12	0.08	102	1.64	0.104
relig_std	0.59	0.13	102	4.50	<0.001	*
scale1	0.03	0.05	102	0.58	0.56
abs:mm	0.00	0.08	102	-0.05	0.961
abs:relig_std	0.12	0.17	102	0.70	0.488
mm:relig_std	0.06	0.13	102	0.42	0.675
abs:scale1	0.00	0.05	102	0.08	0.939
mm:scale1	0.00	0.05	102	0.01	0.993
relig_std:scale1	0.05	0.09	102	0.52	0.604
abs:mm:relig_std	0.03	0.12	102	0.27	0.788
abs:mm:scale1	-0.05	0.05	102	-0.89	0.377
abs:relig_std:scale1	0.00	0.11	102	-0.01	0.996
mm:relig_std:scale1	-0.06	0.09	102	-0.62	0.54
abs:mm:relig_std:scale1	-0.10	0.08	102	-1.23	0.221
China
(Intercept)	0.38	0.13	92	2.86	0.005	*
abs	0.33	0.14	92	2.27	0.026	*
mm	0.45	0.12	92	3.84	<0.001	*
relig_std	0.45	0.10	92	4.50	<0.001	*
scale1	0.03	0.09	92	0.31	0.756
abs:mm	0.09	0.13	92	0.73	0.465
abs:relig_std	0.15	0.10	92	1.40	0.164
mm:relig_std	0.24	0.09	92	2.60	0.011	*
abs:scale1	-0.01	0.09	92	-0.10	0.918
mm:scale1	-0.06	0.08	92	-0.83	0.409
relig_std:scale1	0.00	0.07	92	0.04	0.968
abs:mm:relig_std	0.07	0.09	92	0.71	0.482
abs:mm:scale1	-0.07	0.08	92	-0.86	0.391
abs:relig_std:scale1	0.00	0.07	92	-0.02	0.981
mm:relig_std:scale1	-0.07	0.06	92	-1.27	0.209
abs:mm:relig_std:scale1	-0.04	0.06	92	-0.68	0.498
Vanuatu
(Intercept)	-0.27	0.12	83	-2.32	0.023	*
abs	-0.06	0.15	83	-0.39	0.698
mm	-0.15	0.13	83	-1.13	0.263
relig_std	0.43	0.15	83	2.93	0.004	*
scale1	-0.06	0.09	83	-0.67	0.503
abs:mm	-0.32	0.16	83	-2.00	0.049	*
abs:relig_std	0.19	0.17	83	1.12	0.264
mm:relig_std	0.15	0.14	83	1.02	0.312
abs:scale1	-0.22	0.11	83	-2.03	0.045	*
mm:scale1	-0.03	0.10	83	-0.28	0.783
relig_std:scale1	0.08	0.11	83	0.71	0.479
abs:mm:relig_std	0.28	0.19	83	1.52	0.132
abs:mm:scale1	-0.12	0.12	83	-1.01	0.314
abs:relig_std:scale1	0.25	0.13	83	1.99	0.05
mm:relig_std:scale1	0.02	0.11	83	0.21	0.836
abs:mm:relig_std:scale1	0.18	0.14	83	1.33	0.189

I see two general patterns here:

1. In the US, Thai, and Chinese samples, absorption, porosity, and religiosity all remain significant predictors of spiritual experience when controlling for each other in one of the two analyses reported above. In the other analysis, there is some evidence for “swamping out”: In the US and Chinese samples, the effect of absorption becomes non-significant in the other analysis; and in Thailand, the effect of porosity becomes non-significant in the other analysis. Generally speaking, though, I’d say this is moderate evidence that religiosity does not fully “explain away” the effects of absorption and porosity in these three sites (the US, Thailand, and China).
2. In Ghana and Vanuatu, only religiosity remains a significant predictor of spiritual experience in these analyses. This could be interpreted as evidence that religiosity does “explain away” the effects of absorption and porosity in these two (devoutly Christian) sites (Ghana and Vanuatu).

To help with our interpretation, here are the correlations with religiosity, this time looking at scores standardized within each site:

Joining, by = c("p7_ctry", "scale")

Which experiences are predicted by absorption vs. porosity?

Collapsing across sites

Here are the top ten most strongly correlated items for each of our three predictors of interest (Absorption, Porosity (Dulin), and Martha Story):

Correlation	Question text	Scale
Absorption
0.356	6. I have had the experience of hearing a person's voice and then found that there was no one there. [no one's voice]	Unusual Events
0.326	13. Have you ever experienced a supernatural presence that was not God, a spirit, or a demon? [supernatural presence]	Spiritual Events
0.294	7) [distant msg]. I am completely convinced that it is possible to send a "mental message" to another person, or in some way influence them at a distance, by means other than normal means of communication.	Extrasensory
0.288	1. I hear a voice speaking my thoughts aloud. [voice aloud]	Unusual Events
0.264	12. Have you ever felt a demonic presence as if it was there in the room with you? [Demon in room]	Spiritual Events
0.259	8) [send msg]. I am completely convinced that I have had at least one experience of sending a "mental message" between myself and another person.	Extrasensory
0.258	4. I can hear music when it is not being played. [music]	Unusual Events
0.245	3. I hear people call my name and find that nobody has done so. [call name]	Unusual Events
0.244	5. Have you ever felt that God or a spirit placed an image inside your head? [image in]	Spiritual Events
0.230	19. Have you ever had the experience of being awake but unable to move? [awake paralysis]	Spiritual Events
Porosity (Dulin)
0.623	10. I feel God's love for me, directly. [direct love]	Daily Spiritual Experiences
0.623	9. I feel guided by God in the midst of daily activities. [guided daily]	Daily Spiritual Experiences
0.606	11. I feel God's love for me, through others. [love through others]	Daily Spiritual Experiences
0.575	3. Have you felt that God or a spirit placed thoughts inside your head? [placed thought]	Spiritual Events
0.566	8. I ask for God's help in the midst of daily activities. [God help]	Daily Spiritual Experiences
0.560	4. I find strength in my religion or spirituality. [spiritual strength]	Daily Spiritual Experiences
0.548	1. I feel God's presence. [God presence]	Daily Spiritual Experiences
0.538	12. I am spiritually touched by the beauty of creation. [touched by beauty]	Daily Spiritual Experiences
0.515	3. During worship, or at other times when connecting with God, I feel joy which lifts me out of my daily concerns. [no daily concerns]	Daily Spiritual Experiences
0.500	5. I find comfort in my religion or spirituality. [spiritual comfort]	Daily Spiritual Experiences
Martha Story
0.470	23. Have you ever experienced the presence of God through your own miraculous healing? [heal you]	Spiritual Events
0.465	9. I feel guided by God in the midst of daily activities. [guided daily]	Daily Spiritual Experiences
0.451	10. I feel God's love for me, directly. [direct love]	Daily Spiritual Experiences
0.448	11. I feel God's love for me, through others. [love through others]	Daily Spiritual Experiences
0.443	3. Have you felt that God or a spirit placed thoughts inside your head? [placed thought]	Spiritual Events
0.427	22. Have you ever experienced the presence of God in a miraculous healing (that you saw in person, not on television)? [live heal]	Spiritual Events
0.414	21. Have you ever experienced the presence of God through illness (including as warning or punishment)? [illness]	Spiritual Events
0.413	3. During worship, or at other times when connecting with God, I feel joy which lifts me out of my daily concerns. [no daily concerns]	Daily Spiritual Experiences
0.411	1. I feel God's presence. [God presence]	Daily Spiritual Experiences
0.404	4. I find strength in my religion or spirituality. [spiritual strength]	Daily Spiritual Experiences

And here’s a visualization of the correlations of these three predictors with all “experience” items in this packet (spiritual and secular):

Considering each site individually

Here are the top three most strongly correlated items for each of our three predictors of interest (Absorption, Porosity (Dulin), and Martha Story):

Site	Correlation	Question text	Scale
Absorption
US	0.515	6. I have had the experience of hearing a person's voice and then found that there was no one there. [no one's voice]	Unusual Events
	0.491	13. Have you ever experienced a supernatural presence that was not God, a spirit, or a demon? [supernatural presence]	Spiritual Events
	0.380	3) [psychic]. I am completely convinced that I am psychic.	Extrasensory
Ghana	0.333	1. I hear a voice speaking my thoughts aloud. [voice aloud]	Unusual Events
	0.262	14. I feel a selfless caring for others. [selfless care]	Daily Spiritual Experiences
	0.246	6. I have had the experience of hearing a person's voice and then found that there was no one there. [no one's voice]	Unusual Events
Thailand	0.403	6. I have had the experience of hearing a person's voice and then found that there was no one there. [no one's voice]	Unusual Events
	0.377	1. I hear a voice speaking my thoughts aloud. [voice aloud]	Unusual Events
	0.345	12. I am spiritually touched by the beauty of creation. [touched by beauty]	Daily Spiritual Experiences
China	0.366	15. Have you ever had a feeling of overwhelming emotion during prayer? [emotional prayer]	Spiritual Events
	0.357	14. Have you ever had an experience of uncontrollable shaking or trembling during prayer, or been slain in the spirit? [shaking prayer]	Spiritual Events
	0.352	7. I see shadows and shapes when there is nothing there. [shadows]	Unusual Events
Vanuatu	0.397	13. Have you ever experienced a supernatural presence that was not God, a spirit, or a demon? [supernatural presence]	Spiritual Events
	0.313	3. I hear people call my name and find that nobody has done so. [call name]	Unusual Events
	0.253	5. Have you ever felt that God or a spirit placed an image inside your head? [image in]	Spiritual Events
Porosity (Dulin)
US	0.580	12. Have you ever felt a demonic presence as if it was there in the room with you? [Demon in room]	Spiritual Events
	0.579	1) [esp exists]. I am completely convinced that ESP exists. [ESP is extrasensory perception, or psychic ability]	Extrasensory
	0.554	9. Have you ever had a dream you felt was sent by God or a spirit? [dream]	Spiritual Events
Ghana	0.353	14. Have you ever had an experience of uncontrollable shaking or trembling during prayer, or been slain in the spirit? [shaking prayer]	Spiritual Events
	0.340	9. Have you ever had a dream you felt was sent by God or a spirit? [dream]	Spiritual Events
	0.314	3. During worship, or at other times when connecting with God, I feel joy which lifts me out of my daily concerns. [no daily concerns]	Daily Spiritual Experiences
Thailand	0.479	16. Have you ever had a sense of intense power shoot through you during prayer? [powerful prayer]	Spiritual Events
	0.413	7) [distant msg]. I am completely convinced that it is possible to send a "mental message" to another person, or in some way influence them at a distance, by means other than normal means of communication.	Extrasensory
	0.378	10. I feel God's love for me, directly. [direct love]	Daily Spiritual Experiences
China	0.577	7) [distant msg]. I am completely convinced that it is possible to send a "mental message" to another person, or in some way influence them at a distance, by means other than normal means of communication.	Extrasensory
	0.569	8) [send msg]. I am completely convinced that I have had at least one experience of sending a "mental message" between myself and another person.	Extrasensory
	0.558	1) [esp exists]. I am completely convinced that ESP exists. [ESP is extrasensory perception, or psychic ability]	Extrasensory
Vanuatu	0.444	1. I hear a voice speaking my thoughts aloud. [voice aloud]	Unusual Events
	0.367	12. I am spiritually touched by the beauty of creation. [touched by beauty]	Daily Spiritual Experiences
	0.310	2. I experience a connection to all of life. [connection to life]	Daily Spiritual Experiences
Martha Story
US	0.610	12. Have you ever felt a demonic presence as if it was there in the room with you? [Demon in room]	Spiritual Events
	0.599	13. Have you ever experienced a supernatural presence that was not God, a spirit, or a demon? [supernatural presence]	Spiritual Events
	0.539	7) [distant msg]. I am completely convinced that it is possible to send a "mental message" to another person, or in some way influence them at a distance, by means other than normal means of communication.	Extrasensory
Ghana	0.321	23. Have you ever experienced the presence of God through your own miraculous healing? [heal you]	Spiritual Events
	0.297	7. I see shadows and shapes when there is nothing there. [shadows]	Unusual Events
	0.288	14. Have you ever had an experience of uncontrollable shaking or trembling during prayer, or been slain in the spirit? [shaking prayer]	Spiritual Events
Thailand	0.385	8. Have you ever felt that you tasted God or a spirit? [taste]	Spiritual Events
	0.330	21. Have you ever experienced the presence of God through illness (including as warning or punishment)? [illness]	Spiritual Events
	0.322	16. Have you ever had a sense of intense power shoot through you during prayer? [powerful prayer]	Spiritual Events
China	0.417	13. Have you ever experienced a supernatural presence that was not God, a spirit, or a demon? [supernatural presence]	Spiritual Events
	0.404	1) [esp exists]. I am completely convinced that ESP exists. [ESP is extrasensory perception, or psychic ability]	Extrasensory
	0.389	22. Have you ever experienced the presence of God in a miraculous healing (that you saw in person, not on television)? [live heal]	Spiritual Events
Vanuatu	0.345	1. I hear a voice speaking my thoughts aloud. [voice aloud]	Unusual Events
	0.256	15. I accept others even when they do things I think are wrong. [accept wrong]	Daily Spiritual Experiences
	0.255	2. I experience a connection to all of life. [connection to life]	Daily Spiritual Experiences

And here’s a visualization of the correlations of these three predictors with all “experience” items in this packet (spiritual and secular):

Porosity: good vs. bad emotions

Removed 3 rows containing missing values (geom_point).

Overall scores

Histograms for each scale

Absorption

Daily spiritual experiences

Spiritual events

Ways of being

All items

Excluding reverse-coded items

attributes are not identical across measure variables;
they will be droppedNAs introduced by coercion

Uncommon events

Extrasensory

Hard thinking

All items

Excluding reverse-coded items

attributes are not identical across measure variables;
they will be droppedNAs introduced by coercion

Porosity

Martha story

Demographics

[1] "Participants who passed all attention checks"

[1] "All participants"

[1] "Participants who passed all attention checks"

[1] "All participants"

Important to note that a very substantial majority of participants in Thailand were female, and a slight majority of participants in Vanuatu were female.

[1] "Participants who passed all attention checks"

[1] "All participants"